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Okamoto has obtained a sequence of r-functions for the Pvi system expressed as a double Wron- 
skian determinant based on a solution of the Gauss hypergeometric equation. Starting with integral 
solutions of the Gauss hypergeometric equation, we show that the determinant can be re-expressed 
as multi-dimensional integrals, and these in turn can be identified with averages over the eigenvalue 
probability density function for the Jacobi unitary ensemble (JUE), and the Cauchy unitary ensemble 
(CyUE) (the latter being equivalent to the circular Jacobi unitary ensemble (cJUE)). Hence these 
averages, which depend on four continuous parameters and the discrete parameter A'^, can be char- 
acterised as the solution of the second order second degree equation satisfied by the Hamiltonian 
in the Pvi theory. We show that the Hamiltonian also satisfies an equation related to the discrete 
Pv equation, thus providing an alternative characterisation in terms of a difference equation. In the 
case of the cJUE, the spectrum singularity scaled limit is considered, and the evaluation of a certain 
four parameter average is given in terms of the general Pv transcendent in a form. Applications are 
given to the evaluation of the spacing distribution for the circular unitary ensemble (CUE) and its 
scaled counterpart, giving formulas more succinct than those known previously; to expressions for 
the hard edge gap probability in the scaled Laguerre orthogonal ensemble (LOE) (parameter a a non- 
negative integer) and Laguerre symplectic ensemble (LSE) (parameter a an even non-negative integer) 
as finite dimensional combinatorial integrals over the symplectic and orthogonal groups respectively; 
to the evaluation of the cumulative distribution function for the last passage time in certain models 
of directed percolation; to the r-function evaluation of the largest eigenvalue in the finite LOE and 
LSE with parameter a = 0; and to the characterisation of the diagonal- diagonal spin-spin correlation 
in the two-dimensional Ising model. 
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1 Introduction and summary 

1.1 Setting and objectives 

This paper is the last in a series devoted to a systematic account of the application of the Okamoto T-function 
theory of Painleve equations to the characterisation of certain averages in random matrix theory. The T-function 
theory applies directly to random matrix ensembles defined by a probability density function (PDF) of the form 

1 

-^Ylw{Xj) Yl i^k-Xjf, (1.1) 
j=l l<j<k<N 

where the weight function w{x) is one of the classical forms 



w{x) = < 



e ^ , Gaussian 

x°'e~^ {x > 0), Laguerre ,^ 

a;"(l - a;)'' {0 < x < 1), Jacobi 

{1 + x^y, Cauchy 



The symbol C, which in (1.1) denotes the normalisation, will be used throughout to denote some constant 
(i.e. quantity independent of the primary variables of the equation). The PDFs (1.1) can be realised as the 
joint eigenvalue distribution of Hcrmitian matrices with independent complex Gaussian entries for the first three 
weights of (1.2), and by a stereographic projection of random unitary matrices for the Cauchy weight (with 
rj = N) (see (1.16) below). The underlying matrix distributions giving rise to (1.1) are invariant with respect 
to similarity transformations involving unitary matrices, and for this reason are termed matrix ensembles with 
a unitary symmetry or simply unitary ensembles (not to be confused with unitary matrices). The name of the 
weight function is then prefixed to the term unitary ensemble. In this work our interest is in the Jacobi unitary 
ensemble (JUE) and the Cauchy unitary ensemble (CyUE). We will see that the Cauchy unitary ensemble is 
equivalent to the circular Jacobi unitary ensemble (cJUE) in which the eigenvalues are on the unit circle in the 
complex plane. The averages of interest are 



and 



FM{s;n):={l[{s-Xir) (1.4) 
1=1 



(in the latter, for fi ^ Z a. suitable branch must be specified). 
For AT = 1 (1.3) and (1.4) read 



/S POO 
(s — xYw{x) dx, Fi (s; fi) = (s — xyw{x) dx. (1.5) 
-oo •' —oc 
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A fundamental fact is that as a function of s, and after multiplication by a suitable elementary function of s, these 
functions satisfy a classical second order linear differential equation (Hermite- Weber equation in the Gaussian 
case [29], confluent hyporgeometric equation in the Laguerre case [30], and, as will be shown below. Gauss 
hypergeometric equation in the Jacobi and Cauchy cases). In the Okamoto T-function theory of Piv [53], Py [55] 
and Pvi [54], these same linear differential equations respectively characterise the first member of an infinite 
sequence of r- functions (the zcroth member, as with the random matrix averages (1.3) and (1.4), is unity). The 
general A'^th member of the r function sequence is characterised by the fact that its logarithmic derivative satisfies 
a second order second degree differential equation. Furthermore the A'^th member can be written explicitly as a 
Wronskian determinant and we have shown in the Gaussian case in [29], and in the Laguerre case in [30], that 
the Wronskian determinant is just a rewrite of the average (1.3) (or (1.4) as appropriate). This way we have been 
able to characterise (1.3) and (1.4) in terms of the solution of second order second degree differential equations 
from the Painleve theory. Backlund transformations of quantities associated with the Hamiltonian formalism 
of the particular Painleve system have also allowed us to characterise these quantities in terms of difference 
equations related to discrete Painleve equations. In this work we will show that the same strategy allows (1.3) 
and (1.4) in the Jacobi and Cauchy cases to be characterised similarly. Scaled limits of these averages can then 
be characterised as solutions of limiting forms of the differential equations. Applications are given to the exact 
evaluation of eigenvalue spacing distribution functions, probabilities of last passage times in directed percolation 
models, and the diagonal spin-spin correlation of the two-dimensional Ising model. 



1.2 Definition of averages in the various ensembles 

In the Jacobi case the explicit form of (1.3) and (1.4) is 

Ei,{s;a,b,ix) = ^ dxixlil - xxf^s - xiY ■ ■ ■ l dxN x'}f{l - xn)''{s - xn)'^ {xk-XjY 

_ ^JV(a+/i+JV) 1 / (lxiXi{l — SXlf {1 — XlY ■ ■ ■ i dXN X%{1 — SXN)''il — Xn)'^ 

^ Jo Jo 

n (^fc - (1-6) 



X 

l<j<k<N 



F^{s;a,b, ij,) = ^ j dxi xl{l - xi)'' {s - xiY ■ ■ ■ j dxN x%{l - xn)''{s - xnY i^k-XjY, (1.7) 

^ •^0 •^0 l<j<k<N 

while the explicit form of (1.3) and (1.4) in the Cauchy case is 

E%^s■,v,^^) =^r dx, ii-^ ---r dx^ ii^gii^ n (-^ - -^•)' (i-^) 

O {^+X^)l + ^JV-" l<i<fc<]V 

(a; .) = ^ y_ dx, dx. ^^^Jx, ■ d-^) 

A single quantity which combines both (1.6) and (1.7) is 

Eli{s;a,h,ix;£):=^(^j j ^dxix\{l - xif{s - x\_Y ■ ■ ■ j j ^dxN x%{l - xn)''{s - xn 



r 

X Yl {Xk-Xjf. (1.10) 

l<j<k<N 

Similarly, a single quantity which combines both (1.8) and (1.9) is 

As will be revised in Section 5.1, the power series expansions in ^ of (1.10) and (1.11) are quantities of relevance 
to conditional spacing distributions in the corresponding random matrix ensembles. 
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The integrand in (1.11) requires rj to be real for itself to be real. However, if we write 

(1 + x')-" ^{1 + ixy{l - ix)-", 

where fj denotes the complex conjugate of ri, then the integrand remains real for rj complex, giving a meaningful 
generalisation of the original Cauchy ensemble [14]. Doing this, and setting 7] = rji + ir)2 we can generalise (1.11) 
to read 



i?7(s;(^i,'72),M;C) 7^ / / )dxi 



C\j_^ I (1 + ia;i)''i+'''2(l - ia;i)''i-<''2 

^^-^^)\ _ TT (x^-x^f. (1.12) 

Vj-oc h ' (l + ia;jv)''i+"'2(l-zcc;v)''i-''2 ^^/^{^^^ ' " ^ ' 

It was remarked below (1.4) that for (s — xi)^ to be well defined for ^ Z a definite branch must be specified. 
For s real a natural choice is that {s — xiY is real for a;; < sand {s — xiY = e~''^'^\xi — s^^ for a;; > s. In particular 
this shows 

Ei{s;a,b,fj,;0-=^{J^ )da;ia;?(l-a;i)''|s-a;ir--- -C ^xn x^^il - xn)^\s - xn]" 

X n i^k-Xjf, (1.13) 

l<j<k<N 

where ^ := 1 - (1 - $)e"''''' and similarly for (1.11) and (1.12). 

In the opening paragraph it was noted that the Cauchy ensemble results from a stereographic projection of 
the eigenvalue PDF for random unitary matrices. To be more explicit, following [71], consider the ensemble of 
random unitary matrices specified by the eigenvalue PDF 

^nil + ^'l"^' n l^^-^^l'' zi = e^'>\ eie[-n,n]. (1.14) 

(=1 l<j<k<N 

In [71] this was referred to as the circular Jacobi ensemble with unitary symmetry, and denoted cJUE. When 

u>i — (1.14) is realised by random unitary matrices with the uniform (Haar) measure and is then referred to as 
the circular unitary ensemble (CUE), while (1.14) with = 1 gives the eigenvalue PDF of {N + 1) x {N + 1) 
CUE matrices with all angles 6 measured from any one eigenvalue (taken to be ^ = tt). Making the change of 
variable 

e'* = — X = tan - (1-15) 

1-ix 2 ^ ^ 

(note that 9 = ±7r corresponds to a; — > ±oo) shows 



AT 



i]ii+0;|2'^^ n i^fc-^,fdei...de;v=2^(™i)n— -1— ^ n \xj-xkfdxi---dxN, a-w) 

1 = 1 l<j<k<N 1=1 ^ ' ' l<j<k<N 

thus specifying the relation with the Cauchy unitary ensemble. A generalisation of the cJUE eigenvalue PDF 
(1.14) is the PDF 

1 ^ 

±l[e-'''\l + zi\'-' n l^fc-^^l'' zi = e''\ ei€[-n,7:]. (1.17) 

1 = 1 l<j<k<N 

Under the change of variable (1.15) this transforms into the PDF 



1 ^ 1 



2 

J I 



which corresponds to the generalised Cauchy probability density in (1.12). Consequently, changing variables 
according to (1.15) in 

N 

SS/(0;(a;i,a;2),/x;O := ( 11(1 " ^^(l^..))^"'*' le''^""' " e'^T ) (1.18) 
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and making use of (1.16) and the analogue of (1.13) for shows 

E'j^{4>; {oj^,uJ2), ^l■,C)^ (i^l^^N^/^ ^N'is; (iV + a;i + m/2, 0J2),^i; (1-19) 

The second Jacobi ensemble integral in (1.6) shares the feature of the Cauchy ensemble integral (1 11) of being 
related to an average in the CUE. This is done using the integral identity [24] 

/ dt\t\ ^ ■ ■ ■ I dtN fp^^ f{ti, . . . ,tN) 
Jo Jo 

= (^J^) / dsi e'"""i • • • / da;;ve^"*^"'^/(-e^"*"S...,-e'"'"'') (1.20) 
VsmTre/ J_y2 J-1/2 

valid for / a Laurent polynomial, and making use of Carlson's theorem. One finds 

1=1 " \ 1 I 

where a' = N + a + b,b' = -{N + a) and 

/I/2 /.I/2 N 

dxi--- da;^ n^^^-'^/'ll + n ^(:=e'"*"' (1.22) 

-1/2 J-1/2 i^-^ l<j<k<N 

(note that the left hand side of (1.21) corresponds to the second integral in (1.6) after interchanging /x and b). 
The normalisation (1.22) results from the Jacobi ensemble normalisation 



Jjv(a, 6):= / dxi xj(l — xi)*" • • • / dxN x%{l — xn)*' ]^ {xk — Xj)^ . (1-23) 

•'O •'0 l<j<k<N 



Both (1.22) and (1.23) have gamma function evaluations, 

with the latter following from the Selberg integral [62], and the former by applying (1.20) to the evaluation (1.25). 

The Jacobi average (1.10) and the circular Jacobi average (1.18) (which according to (1.19) is equivalent to 
the Cauchy average (1.11)) admit various scaled N 00 limits. There are three distinct possibilities in that 
the limiting averages can correspond either to the soft edge, hard edge or a spectrum singularity in the bulk. 
The averages at the soft and hard edge have been studied in our previous papers [29, 30], giving rise to Pu and 
Pin transcendents respectively, and will not be discussed here. Instead, attention will bo focussed on the scaling 
to a spectrum singularity in the bulk. This results by replacing X 1— > X/N in (1.18), making a suitable choice 
of the constant C, and taking the N —>■ 00 limit. The problem with C can be avoided by taking the logarithmic 
derivative with respect to X, leading us to consider the scaled quantity 

, N 

uiX; (u;i, u;2), n; {) := ^Imi^ log {Ui^ - T x[lx/^,.))e"^'' |1 + «i|'"Ne'^""''^^^ - ^'l"^)^^^ (1-26) 

(^ is multiplied by X for later convenience). 

1.3 Summciry of the characterisation of the averages as Painleve cr-functions 

Previous studies [65, 35, 72, 12] have characterised Ej^{s;a,b,fi — 0;^) in terms of the solution of nonlinear 
equations related to Pyi . The nonlinear equations obtained have been third order [65, 35] and second order 
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second degree [35, 72, 12]. By the work of Cosgrove and Scoufis [16], these equations are equivalent, and are in 
fact examples of the so called Jimbo-Miwa-Okamoto cr-form of the Pyi differential equation, 

2 2 " 

h' (t{l - t)h"'j + (^/i'[2/i-(2t- l)/i'] + 61626364) ='[[{h' + bi). (1.27) 

k=l 

This fact is most explicit in the work of Borodin and Deift [12] who have shown that 

o-(t) := -t{t - 1)^ log Ej,{l - t; 6, a, 0; i) + bibit + ^(-6162 + 6364) (1.28) 

with 

6i = 62 = iV+^, 63 = ^, 64 = ^ (1.29) 

satisfies (1.27). The quantity (1.11) in the case = has been similarly characterised. In particular, we know 
from [71] that 

a{s) := (l + s')^log^^''(s;(a + iV,0),0;0 (1.30) 

satisfies the equation 

(1 + s^fia'f + 4(1 + s'){a'f - ?.sa{a'f + Acr'^{a' - a) + iasaa' + A[N{N + 2a) - as"] {a'f = 0. (1.31) 
To relate (1.31) to (1.27), change variables 1 1— > (is + l)/2, h(t) t-^ ^h{s) in the latter so it reads 



/i'((l + s'^)h"y + 4(^h'{h - sh') - 161626364)^ + 4 [^(/i' + bl) = 0. (1.32) 

fc=i 

With 

h = (7 — a^s, b = (— o, 0, iV + o, o), (1.33) 



(1.32) reduces to (1.31). 

In this paper we generaUse the Pvi cr-function characterisations (1.28) and (1.30). Consider first the Jacobi 
case. In Proposition 13 we show that 



with 



Uh{t;a,b,^i;(,) := e^[b]t - -e2[b] + t(t - 1)- log S^(t; a, 6, m; (1-34) 

2 at 

b= (i(a + 6) + iV,i(6-a),-i(a + 6),-i(a + 6)-iV-/x) (1.35) 

satisfies (1.27) (the quantities 62 [b] and 62 [b] are defined in Proposition 1). This characterisation can be made 
unique in the cases 5 = 0, 1, when we have the boundary conditions 

lJi,[t;a,b,n;Q) ~ (e2[b] + Ar/i)t + 0(1) (1.36) 

\t\^oo 

U}, (t; a, 6, /i; 1) -^62 [b] -N{a + n + N) + (e^ [b] + iV(iV + o + //) + ^N^^^^^) t (1.37) 
((3.12) and (3.11)). Also, when fj, = with ^ general wc have 

Ui{t; a, 6, 0; - e2[b]t + ie2[b] ~_ C(t - l)p\t) 

o'(t) ^ (1 t)" n^ + b + N+mb + N+l) 

p (t) (1 - 1) r(iv)r(a + iv)r(6 + i)r(6 + 2) ^^■^^> 

where p'^it) denotes the eigenvalue density in the JUE. Regarding the Cauchy case, in Proposition 15 we show 

U^^t; (r?,0),M;O = (t' + 1)| log [{t' + l)^^W/^£°''(t; (r?,0),M;O), (1-39) 

with 

h = {N-r],0,r],-n + i]-N) (1.40) 
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satisfies (1.32). Proposition 15 in fact contains the characterisation of the more general Cauchy average (1.12). 
It follows from this and (1.19) that 

Mjv(wi, wi) 

xexp|-i /V[/°''(cot^;(a;,a;2),M;e)+i(e2[b]-e2[b])-(eUb] + iVA*)cot J)| do] (1.41) 

(eq. (3.25)) where U^^{t;{r)i,r)2),n;^) satisfies (1.32) with 

b = (AT - i]i,ir]2,r]i, -/j. + rji- N). 

The scaled limit (f> >—> X/N, A/" — > oo of the logarithmic derivative of (1.41) is essentially (1.26). This gives 
rise to the differential equation [55] 

4 

{th'{rf - {hv - th'v + 2{h'vff + 4 [h'v + Wfe) = (1.42) 

fc=i 

(c.f. (1.27)) where 

vi+V2+v-i, + Vi = Q (1.43) 

satisfied by a particular auxiliary Hamiltonian hv in the T-function of Py . We show (Proposition 19) that the 
scaled average (1.26) is such that 



satisfies (1.42) with 



h{t) = u{it; (wi, W2),m;S) + + 2wiM + ^^2/2 



ti02 

T 



(1.44) 



2 Overview of the Okamoto r-function theory of Pyi 
2.1 The Jimbo-Miwa-Okamoto cr-form of Pyi 

The sixth Painleve equation Pyi reads 

~ 2\q^ q-t)^'^' t - 1^ q - t)'^ ^ f^{t-lY \^ ^ {q - If ^ {q - tf ) ' ^ ' 

It has been known since the work of Malmquist in the early 1920's [43] that (2.1) can be obtained by eliminating 
p from a Hamiltonian system 

, dH , dH 

In the notation of [41] the required Hamiltonian can be written 

t{t - 1)H = q{q - l)(g - t)p' - [04(9 -l){q-t) + asqiq - t) + (ao - l)q{q - l)]p + a2(ai + a2){q - t), (2.3) 
where the parameters ao, ■ ■ ■ ,04 in (2.3) are inter-related by 

ao + ai + 2a2 + as + 0:4 = 1 (2.4) 
and are related to the parameters a, . . . , 5 in (2.1) by 

a^^ai, /3 = -ia|, 7=^03, 5=i(l-ao)- (2.5) 

One sees that the Hamiltonian can be written as an explicit rational function of the Pyi transcendent and its 
derivative. This follows from the fact that with the substitution (2.3), the first of the Hamilton equations is linear 
in p, so p can be written as a rational function of q, q' and t. The sought form of H then follows by substituting 
this expression for p in (2.3). 

After the addition of a certain linear function in t, the Hamiltonian (2.3) has the crucial feature for random 
matrix applications of satisfying the second order, second degree differential equation (1.27). 
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Proposition 1. [38, 54] Rewrite the parameters ao, . . . , 04 of (2.3) in favour of the parameters 

6i = i(a3 + a4), 62 = i(a4-a3), 63 = ^(ao + ai - 1), 64 = ^(ao - ai - 1), (2.6) 
and introduce the auxiliary Hamiltonian h by 

h = t{t - 1)H + e'2lh]t - ^e2lh] 

= t(t - + (6163 + 6164 + &3&4)t - i Yl (2.7) 

l<j<k<4 

where e'j [h] denotes the jth degree elementary symmetric function in hi , 63 and &4 while ej [b] denotes the jth degree 
elementary symmetric function in 61 , . . . , 64 . The auxiliary Hamiltonian satisfies the Jimbo-Miwa-Okamoto a -form 
ofPvj , (1.27). 

Proof. Following [54], we note from (2.3), (2.2) and (2.6) that 

h' = -q{q - l)p^ + {61 (2g - 1) - 62 }p - 6?, (2.8) 

or equivalently 

q{q - l)ih' + bj) = -(g(g - l)p)' + (bi(2g ~ 1) - b2)q{q - l)p. (2.9) 

We see from (2.8) and (2.9) that a differential equation for h will result if we can express q and q{q— l)p in terms 
of h and its derivatives. For this purpose we note from (2.7) and (2.8) that 

h-th' = q{-h' + e^[b]) - (63 + b4)q{q - l)p - ^e2[h], (2.10) 

while differentiation of this formula and use of the Hamilton equations shows 

t{t - l)h" = 2q{e[[h]h' - e'slh]) - 2q{q - l)p{h' - 6364) - ei[h]h' + eslh]. (2.11) 

The equations (2.10) and (2.11) are linear in q and q{q — l)p. Solving for these quantities and substituting in 
(2.9) gives (1.27). □ 
The T-function is defined in terms of the Hamiltonian by 

//=^logr(t). (2.12) 

In terms of h, it follows from (2.7) that 

h = t{t- 1)^ log ((t - l)"^'''l-^'^=[''lt5«^[''lr(i)). (2.13) 

2.2 Backlund transformations and Toda lattice equation 

Backlund transformations 

T{h;q,p,t,H) = {h;q,p,t,H) 

are birational canonical transformations of the symplectic form. Thus the Hamilton equations are satisfied in the 
variables {h;q,p,t,H). Because there are particular T possessing the property 

TH = h\ , (2.14) 

Backlund transformations allow an infinite family of solutions of the Pyi system to be generated from one seed 
solution. 

Okamoto [54] identified the aiiine Weyl group Wa{D^^^) as being realised by a set of t-invariant Backlund 

transformations of the Pvi system (if Backlund transformations altering t are permitted, a realisation of the affine 
F4 reflection group is obtained [54]). The group Wa{D''^^) is generated by the operators so, . . . ,S4 obeying the 
algebraic relations 

[siSj)""'^ =1, 0<i,j<4 (2.15) 
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where 



[mij] 



(2.16) 



The entries rriij are related to the Dynkin diagram for the afRne root system 





(1) 




Figure 1: D^^ Dynkin diagram. 

Thus for i 7^ i, (0 < i < jf < 4), mij = 2 if the nodes i and j are not connected, while rriij 
connected. An equivalent way to specify the algebra (2.15) is via the relations 



3 if the nodes are 



2 

Si 



1 (0 < «,i < 4), s,sj = sjs, 7^ 2), {siSif = {s2S,f = 1 (i / 2). 



(2.17) 



The operators so, . . . , S4 are associated with afiine vectors qo, . . . ,04 (ao + Qi + 2q2 + Q3 + 04 = 1, with the 
a's regarded as coordinates) in a four dimensional vector space such that Si corresponds to a reflection in the 
subspace perpendicular to a;, and thus Siai = —ai. The action of Si on the other afiine vectors is given by 



where A = (ojj) is the Cartan matrix 



A = 



Si{aj) 



) = 


aj - a 


ttij 







-1 








2 


-1 








I— 1 


2 


-1 


-1 





-1 


2 








-1 





2 



(2.18) 



(2.19) 



(the off diagonal elements in (2.19) are obtained from those in (2.16) by replacing all 2's in the latter by O's, and 
all 3's by — I's). As first identified by Okamoto [54], the operators Si are t-invariant Backlund transformations for 
the Pvi system with action on the parameters specified by (2.18). 

It remains to specify the action of the Si on p and q. For this purpose the most systematic way to proceed is 
to make use of recent work of Noumi and Yamada [50] (see also their subsequent work [51]), who give a symmetric 
formulation of the Backlund transformations for Painleve type systems. In the general formalism of [50], the 
action (2.18) has as its counterpart the action 



Siifj) = fj + 



where U = [uij] is the orientation matrix defined by 






U = 



Cti 



(2.20) 



(2.21) 



The key points are that (2.17) is realised by (2.20) with the fj specified in terms of p, q, t by 

fo = q-t, h=q-oo, /2 = -p, /3 = g - 1, /4 = g. 



(2.22) 
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So 

S3 



-OiO 
ao + Q!2 

ai 

(\:>, 



ai 
ai 

-ai 
ai + a2 

ai 

ai 

ao 



a2 
a2 + ao 

a2 + ai 
— a2 

a2 + as 

a2 + a4 

a2 



as 
as 

as 
as + a2 

-as 

as 

'In 



a4 
a4 

a4 
a4 + a2 

a4 

— a4 

as 



P- 



P 

ao 
q-t 
P 
P 

as 
a4 



pjq - tf + 0:2(9 - t) 
t{t - 1) 

I 1 



q 

q 

ai 
9+ — 
P 

q 

q 

& 

t 



Table 1: Backlund transformations for the Pyi Hamiltonian (2.3). 

(the quantities subtracted from q are the location of the fixed singularities of Pvi ), and that the transformations 
(2.18) and (2.20) together are then Backlund transformations for Pvi . Following [41], the action of the Si on the 
aj , p and q is summarised in Table 1. 

The symmetric formalism of [50, 52] naturally extends the Backlund transformations from a realisation of 
Wa(D^l^^ to a realisation of Wa{D''^'' xi il), where denotes particular diagram automorphisms of d'^^K The 
latter are operators ri and rs (another natural diagram automorphism is r^a = (04, as, 02, ai, ao) and is related 
to r\ and rs by = rirs) with the actions 

ria = (01,00,02,04,03), raa = (03,04,02,00,01) 

of interchanging outer pairs of vertices in the diagram of Figure 1. These operators obey the algebraic relations 

r? = 7-3 = 1, riS2 = S2ri, r3S2 = S2r3 nrj = r^r^, i j, 



Sl = riSoT-i, S4 = riS3ri, S3 = TsSors, 



S4 = r3Sir3. 



Their action on the p and q is given in Table 1. Making use of Table 1 the action of the fundamental operators 
on the Hamiltonian 

t{t - 1)H =: K (2.23) 

as specified by (2.3) is found to be 

soK = K- ao *^*2^^ + ""("3 ^ + ""("4 ~ 1)(^ ~ ^) 
q t 

siK = K 

S2K = K + 02(1 + Ql — ao)t — 02(01 + Q2 + Q3) 

S3K = K — 03(1 — ao)t 
SiK = K- 04(1 - ao){t - 1) 

nK = K — q{q— l)p — 02? + 02(01 — ao)t + 02(00 + 02 + 04) 
rsK = K + {1- t)qp + 02(00 + 02 + a4)(l - t) 

(the first five of these equations can be found in [68]). 
Consider the composite operator 

T-i ■- riSoSiS2S3S4S2, (2.24) 
which from Table 1 has the action on the o parameters 



T301 = (oo + 1, Ol + 1, 02 — 1, 03, 04) 



(2.25) 
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or equivalently, using (2.6), the action on the b parameters 

Tab = (61,62,63 + 1,64) (2.26) 

(in [54] this was referred to as the parallel transformation £3). With K specified by (2.23) we see from (2.24) and 
Table 1 that 

T3K = k\ =K-q{q-l)p-{ai+a2)iq-t), (2.27) 

which was derived in [54] in a different way. There are only three other fundamental shift operators that share 
the property (2.27). They have the actions 

TioK := risiS2S3S4S2SiK = K\ = K — q{q — l)p — a2{q — t), 

IcKl— *(aQ + l,CKl — l,a2,CK3,CK4) 

Tm)K := r3S3S2SiS4S2S3K ~ K\ = K — (t — l)qp, 

I ai— ^(ao+l,ai ,062,0:3 — l,a4) 

T40K := r4S4S2SiS3S2S4K = K\ = K — t{q — l)p, 

I at— »(ao+l,ai ,0(2,02,04 — I) 

although we will not develop the theory of these cases in this work. 

The result (2.27) motivates introducing the sequence of Hamiltonians 

T3H = h\ 

I CXI— »(ao+n,ai+n,a2 —n,cX3, 04) 

and the corresponding sequence of r-functions specified by 

TSH = j^ log r3 [n] , rg [n] = rg [n] (t) = ra {t; 61 , 62 , 63 + n, 64) . (2.28) 

Okamoto [54] proved that T3 [n] satisfies the Toda lattice equation. We will give the derivation of this result using 
the strategy of Kajiwara et al. [41]. 

Proposition 2. The r-function sequence (2.28) satisfies the Toda lattice equation 



r3[n — l]ra[n + 1] 

f3[n] ' " dt 



5'logf3[n]= ''^"-_rj'^'^ , 5 = t{t-l)l (2.29) 



where 

(n+6i+b3)(n+()3+64)/2 



f3[n\ = (t{t-l)) T3[n]. (2.30) 

Proof. Prom the definitions 

5 log r3[n-Jr^[n + l] ^ ^^^^^^^ _ ^^^^^ 

= -q{q - i)p - (ai + 02)(g - 1) 

+ T3-^(qiq-l)p + iai+a2){q-t))\ (2.31) 

where the second equality follows from (2.27). Now T3 is specified by (2.24), and each of the elementary operators 
is an involution so 

= S2S4S3S2SlS0''l = S2S3S4S2SoSiri (2.32) 

where the second equality follows from the commutation relations S3S4 = S4S3, sosi = siso implied by (2.17). 
Using (2.32) and Table 1 we compute that 



-'((?(q-l)p+(ai+a2)(g-t)) 



CK2 (1 — CKO ) 

= q{q - l)p + 0:2(9 - 1) + (1 - cto)q H 1- (ai + «2)t - (cci +02+ 03) 

, / , ^ , „ wi „ n1 g((? - l)p + Q2g + Q2(g - 1) + aj/p „„^ 

+ (a2+Q;3 + Q!4)(l-Q;o) 7 rr -, -T -, ; ; r^- (2.33) 

p q[q - l)p - asq - 04(9 - 1) - a2(a;2 + 03 + Q4)/p 
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Substituting this in (2.31), we can verify that the resuhing expression can be written as 

5 log ^q{q - l)p^ - [osg + 04(9 - l)]p - 02(0:2 + as + a4)j • (2.34) 

To do this we first compute the derivative in terms of p and q using the Hamilton equations, and then check that 
the expression so obtained is indeed equal to the rational function in p and q obtained by substituting (2.33) in 
(2.31) using computer algebra. But according to (2.3), (2.34) is equal to 

<51og [^Js: + Q2(l-ao)]. (2.35) 
Substituting (2.12) in (2.35) and equating the resulting expression with the left hand side of (2.31), we see that 



T3[n - l]T3[n + 1] d r 1 , / ^/1 ^ 

A ~ ■;jT0logr3[nJ + (a2 - n)(l - ao - n), 



where A 7^ is arbitrary. Recalling (2.6) to replace the a's by the 6's and choosing A = 1 gives (2.29). □ 

2.3 Classical solutions 

An identity of Sylvester (see [46]) gives that if 

f3[0] = 1, (2.36) 

then the general solution of (2.29) is given by 



n] = det [5'+''f3[l]] . (2.37) 

L ' 'J ,-.fc=0,l n-l 



As noted by Okamoto [54] and Watanabe [68], the solution (2.36) is permitted by restricting the parameters so 
that 

fei + &3 = or equivalently = (2.38) 

(this corresponds to a chamber wall in the afRne D^'' root system). Thus with 02 = we see that (2.3) permits 
the solution 

p = 0, H = (2.39) 

and thus T3[0] = f3[0] = 1. Furthermore, with this initial condition T3[1] is given by a solution of the Gauss 
hypergeometric equation. 

Proposition 3. [54] Let the parameters a in (2.3) he initially restricted by (2.38), then apply the operator T3 so 
that 

T,H[0] = H[l] = ^ logr3[l](t). (2.40) 
The function T3[l](t) satisfies the Gauss hypergeometric equation 

t{l - t)T^'[l]{t) + (c - (o + 6 + l)t)T:^[l](t) - abT3[l]{t) = (2.41) 

where 

a = —ai =61+64, 6 = ao = 1 + 63 + 64, c = ao + 04 = 1 + 62 + 64. (2.42) 
Proof. It follows from (2.27), (2.23) and (2.38), (2.39) that 

T3K[0] = -ai{q-t) (2.43) 
where q is the solution of the Hamilton equation 



Siq-t) = ^\ - t{t - 1) = ai{q - tf + [(qi + 04)* + (qi + a3){t - 1)] (g - t) - aot{t - 1). (2.44) 

op I "/=o 



Rewriting the left hand side of (2.43) in terms of ''"3[l](t)/T3[l](t) according to (2.40), then substituting the 
resulting expression in (2.44) and simplifying gives (2.41). □ 
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Our interest is in particular integral solutions of (2.41), which being a second order linear equation has 
in general two linearly independent solutions. Consider first the solution analytic at the origin — the Gauss 
hypergeometric function 2F1 — written as its Euler integral representation 

2Fi(a, b; c; t) = ^) "'"'(1 - w)=-'-'(l - du. (2.45) 

We then have 

Tsimt) = T3{t; 61, &2, -bi + 1, 64) = 2Fi(6i + 64, 1 - 61 + 64, 1 + 62 + 64; t). (2.46) 
Integral solutions of (2.41) which in general are not analytic at the origin are given by 

fia,b,c;t)(x r u'-''{l-uY-''-\t-u)-'' du (2.47) 
J p 

where p and q arc any of 0, l,t,;y oo (// = 1) such that the integrand vanishes [36]. Forming from this a linear 
combination with (p, q) = (0, 1) and (p, q) = {t, 1) we deduce from this that T3[l](t) = f^{a, b, c; t) satisfies (2.41), 
where ^ ^ 

f{a,b,c;t)<x(^j^ -ij^ )it"-'=(l -«)=-''-'(< -w)-''dw. (2.48) 

Another case of interest is the linear combination of (2.47) with (p, g) = (—00,00) and (p, g) = (t, 00). After 
suitably deforming these contours in the complex plane we deduce that r3[l](t) = f^^{a,b,c;t) satisfies (2.41), 
where 

f''{a,b,c;—j-)oi\^j j j (1 + m)"""(l - m)""'"' (t - u)"" du. (2.49) 

Let F{a, b, c; t) be any particular solution of (2.41). Then recalling (2.38), (2.42) and (2.30) we see from (2.37) 
that 

f3[n] = det[d'+H''^\t - l)''/^F{a,b,c;t)]j,k=o,...,n-i. 

A useful check on further working is to note that if {T3[n](i)} satisfies the Toda lattice equation (2.29) with 
f3[0] = 1, then {^"/^(t — l)~~'"^^f[n](t)}, 7 arbitrary, is also a solution which is given by the determinant formula 
(2.37). Thus 

f3[n] = t-T"/2(t - l)T"/2 det[6'+H^{t - lfF{a, b, c; t)L-.fc=o,...,n-i (2.50) 

where 

A + B = b, A-B = -f (2.51) 

and the right hand side must be independent of 7. Our task is to substitute the particular solutions (2.45), (2.48) 
and (2.49) in the n x n determinant (2.50) and show that these can be reduced to three distinct n-dimensional 
multiple integrals. 

Firstly we consider (2.50) by taking the solution (2.45). 

Proposition 4. Let F{a,b,c;t) be given by (2.45), and let 5 be given as in (2.29). Then assuming the constraint 
(2.51), 

Un] OC _ ^)i,n/2 + n(n-l)/2 dct[2Fi(a - j, b + t, C; t)]j,k=0....,n-l. (2.52) 

Proof. A fundamental differential-difference relation for the Gauss hypergeometric function is 

t^2Fi(o, 6; c; t) = 6(2^1 (a, b+l;c;t)- 2^1(0, b; c; tfj . 

It follows from this that 

s(t^{t-lf2Fi{a,b;c;t)) 

= (b-A + t{A + B-b)y^{t-l)'^2Fi{a,b;c;t)+bt'^{t-lf+\Fi{a,b+l;c;t). 
In the special case of the constraint (2.51), this reads 

6(t^{t - 1)^2^1 (a, 6; c; tfj = Bt^{t - 1)^2^1(0, 6; c; t) + bt^{t - 1)^+^2^1(0, 6 + 1; c; t). (2.53) 
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Another fundamental differential-difference relation for the Gauss hypergeometric function is 
t{l-t)^2Fi{a,b;c;t) = (c- 0)2^1(0- l,b;c;t) + {a - c + bt)2Fi{a,b;c;t). 

This implies 

6(t^{t-lf2Fi(a,b; c;t)) 

= (^c-a-A + t{A + B- b)^t^{t - 1)^2^1 (a, 6; c; t) + {a- c)f'^{t - 1)^2^1(0 - 1, b; c; t), 

which in the case of the constraint (2.51) reads 

d(t'^{t - 1)^2^1(0, b; c; t)^ ={c-a- A)t^{t - 1)^2^1(0, b; c; t) + {a - c)t^{t - 1)^2^1(0 - 1, b; c; t). (2.54) 

The identity (2.53) can be used to cliniiuatc the operation 5'' from the right hand side of (2.50). Thus 
substituting (2.53) in column k, and subtracting B times column k — 1 (fc = n — 1, n — 2, . . . , 1 in that order) 
shows 

det[6'+H^{t-lf2Fi{a,b;c;t)]j,^=o,...,n-i 

= 6"-Met \sH^{t-lf2Fi{a,b;c;t) 6'+''-h^{t-lf+\Fi{a,b + l;c;t)\ ^^o_..._„_i - 

fe=l',...'n-l 

Repeating this procedure on column fc (fc = n — 1, n — 2, . . . , fe'), for each of fc' = 2, 3, . . . , n — 2 in that order shows 

n-l 

det \6'+''t^(t - lf2Fi(a, b; c; t)l = ff (b)i det \sH^(t - l)^+^Fi(a, 6 + fc; c; t)] 

L J j,fc=0,...,n-l L J j,k=0,...,n-l 

To eliminate (5^ from this expression we substitute (2.54) in row j and subtract {c — a — A) times row j — 1 
(j = n — 1, n — 2, . . . , 1 in that order). This shows 



det [j^t^(t- l)^+^Fl(a,6 + fc;c;^) 



j,k=0,...,n-l 



n-l 

1=1 



t'^(t- 1)^+^^1(0, 6 + fc;c;t) 
(5-'-4'*(t - l)^+'=2i^i (a - 1, 6 + fc; c; t) 



3 = 1, ■ 
k=0,. 



Repeating this procedure on row j, j = n — l,n — 2, . . . ,j' for each of / = 3, 4, . . . , n — 2 in that order gives 
det \5'+''t^{t-lf2Fi{a,b;c;t)\ 

L J j,k — 0,...,n — l 

n-l 

= (_l)"(»-i)/2 T\(^b)j{c-a)jdet\t''{t-lf+''2F,{a-j,b + k;c;t)] 

L i j,k=0,...,n-l 

Removing the factor t^{t — 1)^+* from each column and using (2.51) gives (2.52). □ 
Proposition 5. It follows from (2.52) that 

Tz\n](t) := Tz(t;bi,b2,b3 + n,bi)\ 

1 61+63=0 

OC / d«l---/ dUn\{u^+'"'{l-U,p-'"'-"{l-tU,p-'"' W {Uk-Ujf. (2.55) 
•'O •'° i=l l<3<k<n 

Proof. Prom (2.30) we have 

r3[nl(t) = (t(t-l))-"("-^+'"/^f3[nl(t) 

where use has been made of the equation 61 + 63 = from (2.38) and 6 = 1 — 61 + 64 from (2.42). Substituting 
(2.52) then shows 

Ti[n]{t) OC det[2Fi(a - j, 6 + fc; c; t)],-,fc=o,....„-i. (2.56) 
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Substituting the integral representation (2.45) and recalling (2.42) shows 

T3ln]{t) oc t-"("-i>/"det [J^ dM (1 - u)''^-'''-''-\l - tuy-"'-'" I 

Lj+i(l-«.+i)""'"'(l-i«.+irl • (2.57) 

L J j,k=0,...,n — l 



J j,k=0,...,n—l 
.„.)62-i'3-n(^_j^.)-6i-64 



X det 



The integrand can be symmetrised without changing the value of the integral provided we divide by n!. The 
function of the Mi's outside the determinant is symmetric in the Ui's, so symmetrising the integrand is equivalent 
to symmetrising the determinant. We have 

Sym det - m,+i)""'"'=(1 - tu^+iY] 

L J j,k — 0,. .. ,n—l 

n-1 

= Sym n (1 - tuj+^y det kVi(l - uj+r)"-'-''] 

L J j,k=0,...,n—l 

J=0 

= (Asym n (1 - tuj+,y) det \u^+,{l - mj+i)"-^-'^] 

3=0 

= det (1 — tu-j+i)^ det ti^j.! (1 — m,+i)"~^~* 

L J j,fc=0,...,n-l L ■'^ ^ ' ' J j,fc=0,...,n-l 

= n (1 - Wj+i)""' det [(1 - tuj+xf] det [(-J^i±^) 1 (2.58) 

jj-^ \- \ 3,k=0,...,n-l VW-Uj+lJ \j,k=0,...,n-l 

But the Vandermonde determinant identity gives that for any ai, 02, . . . , a™, 

det[a*+i]j,fc=o,...,n-i = {ak — ttj). 

l<j<k<n 

This allows the determinants in (2.58) to be written as products. Substituting the results for the determinant in 
(2.57) gives (2.55). □ 
Next we will show that by choosing F in (2.50) to equal the solution (2.48), a generalised multiple integral 
representation for fa [n] can be obtained. 

Proposition 6. Let F(a,b,c;t) in (2.50) be given by 

F{a,b,c;t) = J^^ ^u^'^il - uy-^-'^it - u)-" du. (2.59) 

Then 

f3[n]{t) oc **-"/"(*- l)''"/'+"("-^'/'det[F(a- J, 6 + A;,c;t)],-,fc=o,....n-i. (2.60) 

Proof. The proof of Proposition 4 shows that the sufficient conditions for reducing (2.50) to (2.60) is that F 
satisfies the differential-difference relations 

t^F{a, b, c; t) = CoF{a, b + l,c;t)- bF{a, b, c; t) (2.61) 

t{l - t)^F{a, b, c; t) = CiF{a - 1, b, c; t) + (C2 + bt)F{a, b, c; t), (2.62) 

independent of the explicit form of Co, Ci , C2 . To derive the form (2.61) we first change variables w 1— > tw in (2.59) 
so it reads 

F(a, 6, c;t) -^^ )w''-=(l - tM)=-''-'(l - «)"" dw. 

Differentiating this shows 

t^F(a, 6, c; t) = (-c + l)F{a, b, c; t) 
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Writing 

= 1 (2.63) 

we see from this that 

t ^F(o, b, c; t) = -bF{a, b,c;t) - {c - b - l)F{a, b + l,c;t), (2.64) 

thus establishing (2.61). 

To establish the structure (2.62), we first multiply both sides of (2.64) by (1 — t) to get 

t{l - t)^F{a, b, c; t) = btF{a, b, c; t) - bF{a, b,c;t) - {c - b - l)F{a, b + l,c;t) + {c- b - l)tF{a, 6 + 1, c; t). 
But 

tF{a,b+l,c;t) = t-''+^(^J^ )w"-'^-'(^-^)(l-to)'=-''-'(l -«)-"(««. (2.65) 

Using (2.63) and the equally simple manipulation 

(1 - «)-"+' = (1 - w)-" - w(l - w)-" 
shows that the right hand side of (2.65) is equal to 

- (^F{a, b, c; t)+F{a-l, b, c; t)^ + (^F{a, b + l,c;t) + F{a - l,b + l,c;t)^ 

and thus 

t{l - t)^F{a,b,c;t) = btF{a,b,c;t) - (c - l)F{a,b,c;t) - {c- b - l)(^F{a - l,b,c;t) - F(o- 1,6+ l,c; t)^ 
But 

- (c - fe - 1) (^F{a - 1, b, c; t) - F{a-l,b + 1, c; t)) 

(c - 6 - i)r=+^ ( J^^' )u''-''-\i - tuy-''-\i - u)-''+^ du 

+ r^+^a - 1) j^^'' ^u^-^il - tM)^-''-^! - w)-" dw 

= (a - c)F{a - 1, 6, c; t) + {a - l)F{a, b, c; i) 

so we have 

t(l - t)^F{a, b, c; t) = btF{a, b, c; t) + (o - c)F(a, 6, c; t) + (a - c)F(o - 1, 6, c; t), 

in agreement with (2.62). □ 
Following the steps from which (2.55) was deduced from (2.52) allows us to deduce from (2.60) the following 
multiple integral representation. 

Proposition 7. We can rewrite (2.60) and so deduce 

T3[n]{t) = T3{t;bi,b2,b3 + n,b4)\^^ ^'^(^J ~^ ^dui---(^J ^du„ 

71 

X nu7''2-(''=*+")(l-Wi)''"^''+"Hi-Wi)"''"+'"^ n (^k-Ujf. (2.66) 

i— 1 l<j<k<n 

Finally we show that by choosing F in (2.50) to be given by the linear combination (2.49), T3[n] can be written 
as an n-dimensional integral having a form different from the previous cases. 



a+l 
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Proposition 8. Let F{a,b,c;t) in (2.50) he given by 

F{a,h,c;t) = -w)"-''-'(t-w)-"rfw, (2.67) 

■where it is required the integrand be integrable in the neighbourhood ofu = t, ±oo but not necessarily w = 0, 1 (the 
path of integration can be deformed around these points). Then 

f3[n](t) oct''"/'(t- l)''"/'+"("-i^/'det[F(o-i,6 + /c,c;t)],-fc=o,....n-i. (2.68) 

Following the steps from which (2.55) was deduced from (2.52) allows us to deduce from (2.68) the following 
multiple integral representation. 

Proposition 9. We can rewrite (2.68) and so deduce 

j jdul--- (^J -i j j dUn 

X f[Mr''^-(''3+»)(l_„.)i'2-(63+n)(^_^.)-(6i+64) ^ {Uk - Ujf . (2.69) 
i=l l<j<k<n 

The only necessary detail of the contours (intervals) of integration in (2.69) is that the integrand vanishes at 
the cndpoints. Deforming the contours in this manner we see from (2.69) that 



T3{ — - — ;6i,62,63 +n, 64) °^ ( / ~W )dui---( / / )dw„ 

Z 161+63=0 ^J-oc Jt ' Jt ' 



X [](l + mi)-''=^-(''3+")(l-mi)''""''+"Hi-«i)''"'' n {Uk-Uif. (2.70) 

i=l l<j<fc<n 

2.4 Schlesinger Transformations 

In this part we develop difference equations arising from the sequence generated by the action of T3, which are 
also known as Schlesinger transformations because the formal monodromy exponents are shifted by integers. In 
doing so we will demonstrate that this recurrence in the canonical variables can be expressed in a form which is 
precisely that of the discrete fifth Painlcvc equation dPv • We thus establish directly that the discrete dPv is the 
contiguity relation of the continuous Pyi equation in contrast to other treatments [61, 60, 33, 45]. The Tz.,T^^ 
operators are equivalent to the J?(9), J?(io) operators respectively in [45]. 

Proposition 10. The sequence {g[n],p[n]}5^Lo generated by the shift operator Ti with parameters a = (ao + 
n, ai +n,a2 — 71,03,0:4) defines an auxiliary sequence, equivalent to the sequence {fflfi], /[njj^o satisfying the 
discrete fifth Painleve equation dPy 

t-1 /[w](/[nj +03) 

f[n]+f[n-l] = -a3 + ^^ + -^^ ^ (2.72) 

gin] - 1 t{g[n] - 1) - g[n] 

where 

9-=—^, f ■= g(g - 1)P+ (1 - 02 - a4)(g - 1) - asq - qo '^^'^_ . (2.73) 

51 q t 

Proof. To begin with we construct the forward and backward shifts in the canonical variables under the action 
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of T-i using Table 1 {q :— g[n],p := p[n]) 

q[n+l] :=q = -L- [{q-l){q-t)p+{ai + a2){q-t)-aoit-l)] 
q I 

X [q{q - - t)p + (ai + 02)9(9 -t)- ao{t - l)q + a^iq - t)] (2.74) 



p[n + 1] := p 



{q ~ t) \{q - t)p + Qi + 02] Xu 

tit - 1) [9(9 - t)p + (ai + 02) (g - t) - aot] [{q - 1) (g - t)p + (qi + 02) (g - i) - ao(i - 1)] 

(2.75) 

q[n -l]:^q= , " 1)P + [9(g " + + ^2.76) 



pr„_l, .__p_ f [{q-t)p + a2] ^ ap - 1 

t{t — 1) \ [qp + 02] [{q — l)p + 02] q{q — l)p^ — (0:4(9 — 1) + 039)^ — 02(02 + 03 



+ 04) } ' 



(2.77) 



where 



Xu = q{q - l)(g - t)^p^ - [04(9 - l){q - t) + 039(9 - t) + {ao - 1 - 01)9(9 - 1)](9 - t)P 

+ (01 + 02)^(9 — t)^ + (01 + 02)[(0l + 02 + 04)* + (Ol + 02 + 03)(t — 1)](9 — t) 

-00(01 + l)t(t- 1) (2.78) 

Xd = q{q -l){q- t)p^ - [04(9 - l){q - t) + a^qiq - t) + (oo - 1 + 01)9(9 - 1)]P 

+ 02(9 -t) + 02(02 + 04)t + 02(02 + 03)(t - 1), (2.79) 

(note that the factors Xd, Xu are distinguished by being quadratic in p). Prom these variables certain products 
can be constructed that have simple factorisable forms 

q P = [(g ~ + 0:1+ Q2] [9(9 - l)(g - t)p + (01 + 02)9(9 -t)- ao{t - 1)9 + 0:4(9 - t)] ,2 80) 

t-1 [q{q - t)p + {ai + a2){q - t) - aot] ^' ' 

(-^ 1 1" ^ [(9 - t)P + 0:2] [9(9 - 1)(9 - t)p + (oi + Q:2)(g - 1)(9 - ^) - aot{q - 1) - 03(9 - t)] 

t [{q-l){q-t)p+{ai+a2){q-t)-ao{t-l)\ ^ ' 



(9 - t)p 



[(9 - t)p + ai + 02] 



[9(9 - *)P + (ai + a2)(9 - i) - oot] [(9 - 1)(9 - <)P + (ai + o:2)(9 - t) - ao(t - 1)] 
X -{ 9(9 - 1)(9 - <) V - [0:4(9 - 1)(9 - i) + 039(9 -t) + 2009(9 - 1)](9 - t)p 



+ ao[o4(9 - 1)(9 - *) + 0:39(9 - *) + 009(9 - 1)] - (1 - O2)(ao + 01 + 02)(9 - tf } (2.82) 



9 - 1 _ 1 9^(9 - l)p^ + [2029(9 - f ) - Q3g]p + ojg - 02(02 + oa) 
^-t t 9(9 - l)p2 - [04(9 - 1) + 039b - 02(02 + 03 + 04) 



(2.83) 



Prom the ratio of (2.80,2.81) one notices that it can be simply expressed in terms of the single quantity / defined 
by (2.73) which is the first of the coupled recurrences for dPy (2.71). To find the other member of the pair we 
evaluate / := f[n — 1] using the above formulas and find 

7 = -(9-l)(9P + 02), (2.84) 

so that when this is added to / we arrive at (2.72). □ 

Although of theoretical interest, Proposition 10 does not immediately lead to a difference equation for the 
Hamiltonian itself. However such an equation can be derived by using the workings from the derivation of 
Proposition 10. 

Proposition 11. The sequence of Hamiltomans {K[n]}^=o generated by the shift operator T3 with parameters 



18 



a = (ao + n,ai + n, 02 — n, 03,04) satisfies the third order difference equation 

1^(1 — ao)K + aoK^ + ai)K — aiK + (ai + a2)(l + ai — ao)t — (ai + a2)(ai +02+ 03)^ 



j^/sT - /sT + 1 - ao - 0:4 - (1 + ai - ao + 2a2)tj 



K — K — ao — 04 + (2ao + 03 + a4)t 



+ t 



(t — — ao — ai)if — ai(l — ao)(-R' — + — ao)(ai+a2)(l — ao + ai)f — (1 — ao)(ai+a2)(ai + a2 + a3)j 

X 1^(1 + ao + ai)K — ao(l + ai){K — K) + ao{ai + a2)(l — ao + ai)t — ao{ai + a2)(ai + a2 + 03 )j = 

(2.85) 



where K := t{t-l)H[n],K := t{t - l)H[n - 1], K := t{t - l)H[n+l],K := t{t - l)H[n + 2]. 

Proof. We focus our attention on the quantity 

Z := -q{q - \)p - (ai + az) (9 - , (2.86) 

which from (2.27) \s Z = K[n + 1] — K[n]. First we consider Z := Z[n — 1]. Prom the definition Z is related to / 
so employing (2.84) and (2.83) for the downshifted variables we have 

Z = {q- l){qp + 02) + (1 - ao)q - 03 + (ai + 02){t - 1) 

(2a2 + as + a4)[g(g - l)p + 029] - 02(02 + 03) ,^ 
q{q - l)p2 - [04 (g - 1) + 03g]p - 02(02 + 03 + 04) ' 

We transform this expression by replacing the canonical variables q,q~ l,p where possible by K, Z and retaining 
q — t. The resulting equation is then one which is a linear equation for q — t m terms of K, Z, Z , 

q-t=[K+{l- ao)Z] + Z + 02 + 03 - (oi + a2){t - 1) - (1 + 02 - oo)t] 

/ [(1 - oo - ai)K - ai(l - ao)(^ + Z) 
+ (1 - ao)[-oi(o2 + oa) + (o? + 01O2 + 02 + a2a3)(t - 1) + (oi + 01O2 - ooOi + al + a204)t]j . (2.88) 

We use a similar strategy and evaluate the upshifted Z := Z[n + 1] using (2.81,2.80,2.82,2.75). We find after 
considerable simplification 



Z = q{q— l)p + (01 + a2){q — t) — (2ao + 03 + 04)* — (1 + 00 + oi) 



q-t 



+ (1 + oi) /*-^ { (1 + ao + ai)q{q - l){q - t)p + (ai + a2)(l + oo + oi)q^ 



{q - t)Xn 

+ [-(1 + 00 + ai)(ao + 01 + 02)* + 00(00 + 04) - (01 + a2)(ai + 02 + 03)]? 

+ (oo+oi+a2)(l-a2-04)t^, (2.89) 

and after the variable replacements described above one has a linear equation for 
t{t-l) 



q-t 



= — ^K — aiZ + (oi + 02)[(1 — Oo + oi)t — oi — 02 — 03]] ^Z + Z — ao — 04 + (2ao + 03 + a4)tj 

(1 + Oo + ai)K + (1 + Qi — aoo\)Z — oo(l + ai)Z 



+ oo(oi + a2)(l — Oo + a\)t — oo(oi + 02)(ai + 02 + 03) > (2.90) 
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By eliminating q — t between (2.88) and (2.90) and further simplification one arrives at the stated result (2.85). 

□ 

Wo remark that the right-hand side of (2.86) can be written in terms of / and g, and then coupled with the 
dPv equations (2.71,2.72) to provide an alternative scheme to calculate K[n]. A similar recurrence appears in 
the recent work of Borodin [10]. 

3 Application to the finite JUE and CyUE 

We are now in a position to identify the multiple integral representations for the r-functions of the classical 
solutions to the Pvi system, (2.55), (2.66) and (2.70), with the spectral averages defined by (1.6), (1.10) and 
(1.12) respectively. 

3.1 The JUE 

The T-function solution (2.46) is relevant to the average (1.6), for with N = 1 we have 

Ei{t; a, b, n) oc t"+^+'2i^i(-6, a + 1; fj, + a + 2;t) 

= e+''+\3{t\ -(a + 6)/2, M + 1 + (o + 6)/2, 1 + (a + 6)/2, (a - 6)/2). (3.1) 

The multiple integral representation (2.55) of T3[n](i) is of the type occuring in the definition (1.6) of -Eiv(t; a, b; fx). 
This allows the latter to be identified with a T-function for the Pyi system, and its logarithmic derivative identified 
with an auxiliary Hamiltonian (2.7) and so characterised as a solution of the Jimbo-MiwarOkamoto cr-form of the 
Pvi equation (1.27). 

Proposition 12. Let Tz[n]{t) = T3{t; 6i, 62, 63-l-n, 64) refer to the T-function sequence (2.28) with rsfO] = 1 

and r3[l](t) given by (2.4.6). Then we have 

Ei{t; a, b, 11) = Ct^^'+^^^<-^'+^''^T3{t; b), b = (-(a + b)/2, fj, + N + {a + b)/2, N + {a + b)/2, (a - b)/2) (3.2) 

and consequently 

tit - 1)^ log ((t - l)'=2[bl-ie2lb]^le2[b]^-(bi+fe3)(62+04)^J^(^. J,^ ^ [;J J,^ ^) (3 3) 

where U^{t;a,b,fji,) satisfies the Jimbo-Miwa-Okamoto a-form of the Pyi equation (1.27) with b specified as in 
(3.2), and e2[b], e2[b] defined as in Proposition 1. The latter is to be solved subject to the boundary condition 

UUf. a, b, ,) ~„ 4e2[b] + (4[b] + m-^±^)t 

4^' 4(M + a - ;')iV + \bia + b) - \Ka - b) + ^-bN -'-ia + bf + ^^^^^tW) ' 

(3.4) 

with Ul) given by a power series in t about t = 0. 

Proof. The first equation follows immediately upon comparing (2.55) with (1.6), and rewriting the exponent 
in the first factor in (1.6) in terms of the 6's. The equation (3.3) then follows from (2.13). For the boundary 
condition, we see from the second integral in (1.6) and the definition of the 6's that 

r^(«+''+^^S^(t; a, b, M) = t-^'^+'^^^'^+'^^EUt; a, b, m) ~ ^^^^ (l - b^^^^^^Pf^^^t) (3.5) 

where J^Ar(a, /u)[Xi^i a:^] is the integral (1.23) with an extra factor of Xi^i iii the integrand. The ratio of 
integrals in (3.5) can be evaluated using a generalisation of the Selberg integral due to Aomoto [3] to give (3.4). 
Alternatively, the 0(1) term in (3.4) can be substituted in (1.27) to deduce the 0(t) term. □ 
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Let us now reconcile Proposition 12 with the result that (1.28) satisfies (1.27) with parameters (1.29). For 
this purpose we note from (3.2) that with = 0, 62 = 63. Then the exponents in (3.3) simplify, and we see that 

Ui,{t; a, h, 0) = -thi + i(6i64 + bl) + t{t - 1)^ log Ej,{t; a, b, 0). (3.6) 

Interchanging a and 6 as required by (1.28) changes the sign of hi, so with b given as in (3.2) with /i = 0, we see 
from (3.6) that 

-U^{1 - t; b, a, 0) = -tbl - i(6i64 - bl) + t{t - 1)^ log £^(1 - t; a, b, 0). 

But in general, if h(t) satisfies (1.27) with parameters b, then —h{l — t) also satisfies (1.27) but with an odd 
number of 61, . . . , 64 reversed in sign. Choosing to reverse the sign of 64, it follows that 

-tbl + ^(6164 + bl) + t{t - 1) log£;^(l - t; a, b, 0) 

satisfies (1.27) with b given by (a suitable permutation of) (1.29), which is in agreement with (1.28). 

A corollary of Proposition 12, which follows from (1.21) and the second equality in (1.6), is the formula 

, N 
UUt;a,l^, b) = t{t - 1)^ log Ut - l)4[b*]-ie.[b*]^le.[b*] / TT ^ ^ ^^^^A 

€tZ \ i 1 CUEjy 

b* = ( - (o + m)/2, b + N+{a + n)/2, N+{a + n)/2, (o - n)/2j . (3.7) 
Comparison of the r-function solution not analytic at the origin (2.48) with (1.10) shows, for A?^ = 1 
Ei{t; a,b,ii;^) oc /""(-/U, -/U - a - 1 - 6, -// + o; t) 

= T3(t; 1 + i(a + b), i(& - a), -i(a + b),-fi-l- i(a + fe)). (3.8) 

Furthermore comparing the n-dimensional integral (2.66) with (1.10) allows us to characterise the latter in terms 
of a solution of (1.27). 

Proposition 13. LetT3[n]{t) = r3(t; 61, &2, fea+^i, ^4) refer to the t -function sequence (2.28) withrzlQ] = 1 

1 61+63=0 

and r3[l](f) given by (2.59). Then we have 

Ej,{t;a,b,n;0 = CT3{t;h), b = (i(a + 6) + A, i(& - a), -i(a + 6), -i(a + &) - iV - (3.9) 
and consequently 

e^[h]t - ^e2[h] + t{t -1)^ log Eir(t: a, b,^i;^)^U'r,{t; a, b,fi;0 (3.10) 

where tj}f{t;a,b,fj,;£^) satisfixis the Jimbo-Miwa-Okamoto a-form of the Pyi equation (1.27) with b specified as in 
(3.9), and e2[b], e2[b] defined as in Proposition 1. In the case = 1 we have the boundary condition 

UUt; a, b, m; 1) -U^ [b] - A(a + m + A) + (ei[b] + {N + ix + a)N + bN^±^—) t 

- ^(M + « - ^)iV + + b) - - b) + {-bN - i(a + bf + bN^f±f^) t 

(3.11) 

with U}f a power series in t about t = 0, while in the case ^ = 

[/^(t;a,6,M;0) ~ ( e'ilh] + N^i)t + 0{1) 

|t|— >oo \ / 

. , y (3.12) 

with U}f, apart from the leading term, a power series in 1/t. 
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Proof. The only remaining task is the specification of the boundary conditions (3.11) and (3.12). The first of 
these follows from the fact that Ej^{t;a,b, fi;!) = E^{t;a,b, fi) and then substituting (3.5) in (3.10), while the 
second follows from the fact that EN{t;a,h, ix\Q) = F^{t;a,b, fi) and noting from (1.7) that F^{t;a,h, n) ~ t^** 
as t — > oo. □ 
Comparing the definition (1.10) with (1.6) shows 

Ei,{t;a,b,^i;l) ^ ENit;a,b,fi) (3.13) 

so the cr-function U in (3.10) must be related to the cr-function U in (3.3). To explore this point, writing b in 
(3.2) as b = (6i, 62, 63, 64) we see from (3.9) that 

b = (63,-64,61,-62). (3.14) 

The differential equation (1.27) is unchanged by the replacement of b by (3.14), so U and U in fact satisfy the 
same equation. In fact U and U are in this case the same function, as the identity (3.13) together with the readily 
verified formula 

e2[b]i - ^e2[b] - {t - l)(6i + 63)(62 + 64) = e^[b]t - ^62 [b] 

show that the left hand sides of (3.3) and (3.10) agree. 

Another special case of (1.10) of interest is the coefficient of 

[^'^]E'i,{s;a,b,fj,;0 = ^-!— dxi x1(l - xi)\s - xiY ■ ■ ■ cLxn xl^il - xn)\s - xnY W {xk-Xjf. 

■'^ l<j<k<N 

This exhibits the functional property 

K^]^^(i;a,6,M;0 = (-l)^s-^("+''+''+~'K^]S^(5;a,6,M;0- 
The r-function evaluation (3.9) then implies 



where 



r3(i;b) = (-l)^t'^3+^*«^i+^^V3(t;b) (3.15) 

b = (i(6i - 62 + 63 - 64), ^(-61 + 62 + 63 - 64), i(6i + 62 + 63 + 64), i(-6i - 62 + 63 + 64)) (3.16) 

(b is obtained from b by simply interchanging 6 and n in the latter; recall (3.14) and (3.2)). This result can be 
understood within the context of the general theory of the Pyi equation. Thus the mapping 

iq,p, H,t;h)^{-, (61 + 63)9 - q^p, - § + $(t), i; b), $(<) = -^(63 + 61) (63 + 64) 

q t-' t t 

has been identified in [54] as a canonical transformation of the Pyi system (the value of $(i) was not given 
explicitly in [54]). Recalling (2.12), up to a proportionality constant this immediately implies (3.15). 

A consequence of the difference equation (2.85) is that a difference equation for f7^(t;o, 6, /x;^) with respect 
to can be found. 

Proposition 14. The Jimbo-Miwa-Okamoto cr-function U}f{t;a,b, iJ,;^,) satisfies a third order difference equation 
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in n, for all N, a, b, t 

\^^{N + a + n + l)U - {N + a + ^i)U + 1/4(0 - bft - 1/4(0'' + 6^) - 1/4(0 + b)N^ 

x[{N + b + fj,+ l)U-{N + b + n)U + 1/4(0 - bft + y^ab + 1/4(0 + 6)iv] 



[f/ - f/ + (2iV + 2/i + o + b)it - 1/2)] 



U-U+{2N + 2n + 2 + a + b){t - 1/2) 



= t{t - 1) [(AT + o + ij){N + 6 + /x)([/ - f/) - (2A/- + 2^ + o + b)U 
- 1/4(0 - 6)^(2Ar + 2^i + o + 6)t + 1/4(6 - a)(2iV6 + (6 - o)/i + 6(0 + 6))] 

X [(TV + o + M + l)(iV + 6 + /i + l)(t? - J7) - (2iV + 2^i + 2 + o + 6)C/ 

- 1/4(0 - 6)^(2Ar + 2(0 + 2 + o + 6)t + 1/4(6 - o)(2iV6 + (6 - o)(/i + 1) + 6(0 + 6))] , (3.17) 

where U ■- U}f{t;a,b, ^■,(,),U := U}f{t;a,b,fj,- l;C),f/ := Ulfit; a,b, fi + 1;^) etc. 

Proof. This follows from (2.85) using the parameters in (3.2) with 62 <-> 63 and the invariance of ft or ?7 under 
this interchange, and the relation 

U},it; o, b, ti; = K[ti] - [b{N + At) + 1/4(0 + 6)'] t - VaiV" - i/2iV(y:i + o - 6) - 1/4^(0 - 6) + 1/46(0 + 6), (3.18) 

and the shifted variants. □ 



One can easily verify that (3.17) is invariant under the transformations o<->6, ti— >1 — t and U}/ 1— > —Ul 



3.2 The CyUE and cJUE 

Comparison of the r-function solution (2.49) with the average (1.12) shows, for A?^ = 1 

E?^it; ('?i,^2),m;0 oc /*^''(-/i,2r7i -/i - 1,771 +ir]2 - m; 
'1 + it 



/1 + It \ 

T3[^—^;-m + 1,«»?2,7?1,7?1 -/i - Ij, (3.19) 



where in obtaining the second line use has been made of (2.42) with 63 = —61. For n > 1 comparison of (2.70) 
with (1.12) allows us to deduce the analogue of Proposition 12 for . 

Proposition 15. Let T3[n]{t) = T3{t; 61, 62, 63 +n, 64) refer to the T-function sequence (2.28) with r3[0] = 1 

lbl+i)3=0 

and r3[l](t) given by (2.67). Then we have 

E%^t;{r,i,r]2),li;0 = CT3C^^;h), b = (AT - 771, 1772, -M + r?i - iV) (3.20) 

and consequently 

(f + i)^iog ((it - l)«^W-^^W(rf+ l)i^-l''l£^^(t;(,7i,,72),Ai;C)) = t^^^(^;(»?i,»72),M;0 (3.21) 

where U'^^ {t\ (j?i, 772), A*; satisfies the transformed Jimbo-Miwa-Okamoto a-form of the Pyi equation (1.32) with 
b specified as in (3.20), and eifb], 62 [b] defined as in Proposition 1. In the case ^ = 1 we have the boundary 
condition 

U^^t;{vi,V2),l.i;l) ~ (e'2[h]+N{N + ^l-2r^l))t 

t— ►-00 \ / 

. _,2,_.2(iV- 770(A^ + M-m)^O(,/,) ^'-''^ 

t->-oo 771 

while tn the case 5 = 

U^''{t;{m,m),l^;0) ~ (e'2[h] + N^i)t + 0{1) 



t — ► — 00 



(3.23) 

{N - r]i)H + 0{1) 
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We see from (3.21) that with 772 = we have 62 = 0. This in turn implies that 62 [b] = 62 [b] and so (3.21) 
reduces to 

it' + 1)| log [{t' + ir'-^^^/'E^^t; (r?i,0),/x;O) = U^' {t; (r?i,0),M;O- (3-24) 
That this with /j. — 0, rji — N + a, h = (—a, 0,N + a, a) satisfies (1.32) is in precise agreement with the fact noted 
earlier that (1.32) with the substitution (1.33) gives the equation (1.31) satisfied by (1.30). 

According to (1.19), we can write down from Proposition 15 an analogous result for {ip; {uji,uj2), fi;^)- 
Using the fact that for — > this quantity tends to a constant, and evaluating the latter in terms of Mjv(o, 6) as 
specified by (1.22), we have 

EN{(t>;iu>i,uj2),n;e. ) = 5^ ^ 

X exp I - i r (U^^cot ^; {N + ioi + /j/2, L02), n; + 0J2{N + wi - 1/2^) + (wi + 1/2/^)^ cot ^dO^. (3.25) 
J 

In the case a;2 = the average E"^ satisfies a functional relation which according to (3.25) must be related to a 

corresponding functional relation of U'^ . Thus, in the case a;2 = the integrand of the multi-dimonsional integral 
specifying E'j^ is periodic, which in turn allows the change of variable 9i 61 — (j). From the latter it follows 

SS/(0;(a;i,O),2/i;O = ^^'(0;(M,O),2a;i;O. (3.26) 

To exhibit this symmetry in (3.25), we note from (3.20) that in the special case a;2 = we must substitute in 
(3.25) 

b = (-(a;i+/[//2),0,iV + a;i+///2,a;i-/x/2), eafb] = e2[b], e2[b] + iV/i = -(wi + ///2)^ (3.27) 

Replacing /it by 2// as required by the left hand side of (3.26) we see that all quantities in (3.27) are then symmetric 
in /i and uji as required by the right hand side of (3.26), except for the component = oj\ — /i/2, which is 
antisymmetric under these operations. However we see the equation (1.27) in the case 62 = is unchanged by 
the mapping 64 1— > —64 and so indeed (3.25) is consistent with (3.26). 

With ^* specified as in (1.13), the average defining .E^(0; (a;i, a;2), 2/i; ^*) is proportional to the CUE average 

1 , 



;na-rx[l,,.))e-*ii+.<r-(^)''(i+*.o^^),^J^^^,,=:^.(t;-i,-2,M;r (3.28) 

It then follows from (3.25) that 

- log AN{t\ wi, W2, m; D 

= \ (t^^^(i^; (Af + cui + ^l,u:2), 2m; + i(ei[b] - 62 [b]) - i^(e^[b] + 2N^i)) (3.29) 
where U'^ satisfies (1.32) with 

b = (-(wi + m), iw2, AT + + - m) =: b'^''- (3.30) 

But, in the case ^* = 0, the same average (3.28) results from (3.7) upon introducing a factor of t~>^^l'^ into the 
average and making the replacements 

6 h- > 2a;i, /x i— » 2/i, a i— » — (iV + + ia»2 + m)- 

The quantity U]^ in (3.7) then satisfies (1.27) with 

b= (^i(Ar + u;-/i),w+i(Ar + w + /;i),i(Ar-w + /i),-/i-i(Ar + w + M)j, w = wi + itJ2. (3.31) 



It follows that 



t (t - 1) — log ^jv (t ; wi , a;2 , m; 0) = f/i, (t ; - (iV + wi + ia;2 + m) , 2m, 2a;i ) - Cit + C2 (3.32) 
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(3.36) 



where 

Ci = 62 [b] + nN, C2 = ie2 [b] + nN. (3.33) 

Comparing (3.29) and (3.32) shows 

([/^"(i^; (Af + c^i + M, c^2), 2m; 0) + i(e^[b] - 62 [b]) - 7^(e'2[b] + 2Nf,)) 

= Ui,{t; -(A/" + tJi + iLJ2 + m), 2m, 2wi) - Cit + C2. (3.34) 
In fact (3.34) is a special case of the following transformation property of (1.27). 
Proposition 16. Let h satisfy (1-27), and put 

hit) = i{t - l)l/(i^) -l(bl+bl+bl + bi)t. (3.35) 

Then f{s) satisfies (1.32) with b 1— » b, where 

hi = 5(61 - &2 + 63 + 64), 62 = ^(61 - 62 - 63 - 64) 
63 = 5(61 + &2 + 63 - 64), 64 = i(6i + 62 - 63 + bi) 

(or any permutation of these values). 

Proof. Substituting (3.35) in (1.27), changing variables to 

.t + l 
s = I 

t- 1 

and equating like terms with (1.32) modified so that h i—^ f and b i—>b, shows the statement of the proposition is 
correct provided 

1 22I I2 

616263&4 = Ye^^^I*' ~ -jbib2h3hA - ;je2[b ] 

-e2[b"] = -l{ei\h']f - ibibibM + ie2[b'] 

O Z 

-e3[b^] = -iei[b'l6i 626364 - ^(ei[b'])' + iei[b''le2[b'] - e3[b']. (3.37) 

Direct substitution of (3.36) into these equations verifies their validity. □ 
We can check immediately that b^^ in (3.30) is related to b £is in (3.31) according to (3.36), and furthermore 
that 

^f^{e'2[h°^] - 62 [b^^]) + ^(1 + i)(e2[b^^] + 2iV/i) + Cit - C2 + ^(6? + 6i + 6^ + bl)t = 0, 
which together show (3.34) follows from Proposition 16. 

An analogous result of the difference equation in /i for the Jacobi case (3.17) holds also for the Cauchy 
Jimbo-Miwa-Okamoto cr-function U'^ {t; (r/i, 772), pi; C)- 

Proposition 17. The Jimbo-Miwa-Okamoto a-function f/^'*'(t; (??i, 772), m; satisfies a third order difference 
equation in n, for all AT, 771, 772, t 



[{N + n+l- 2?7i)[/ -{N + fi- 2ni)U + nit + jfAN - m)] 

X [(AT + M + l)U -{N + n)U + vh -r]2{N- 



X ^U-U + 2{N + fj,-ni)t\ 



U-U + 2{N + n + l-rii)t 



= (1 + t^) [(A- + n){N + M - 27]i){U -U)-2{N + n- m)U -2r,i{N + n- ni)t - 2niri2{N - 771) 
X [(AT + M + l)(Ar + /X + 1 - 277i)(C? -U)-2{N + H + r]i)U - 2rn{N + + 1 - 7?i)t - 27?i772(Ar - 

where U := U^^t;{vi,V2), l^;0,U := U^^t-,{r,i, 7,2), (.i - 1;0,U := U^^t, {7,1,7,2), + etc. 



(3.38) 
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Proof. This follows from (2.85) using the parameters in (3.20) with 63 64 and the invariance of ft. or f7 under 
this interchange, and the relation 

U^^t; {m,r]2), m; = K[fi] - [{N - + HN] t + 1/2 [(AT -r,if + ^iN + irnil^ - 771)] , (3.39) 

and the shifted variants. Finally one has to transform the variables 1 1— » (1 + it)/2 and U'^^ ^ iU'^^ /2. □ 

3.3 Duality relations 

In our previous studies of the average (1.4) in the Gaussian and Laguerre cases we have exhibited a duality in // 
and A'" which in fact extends to random matrix ensembles in which the exponent 2 in the product of differences 
(1.1) is replaced by a continuous parameter (3. The same holds true of the average (1.21). Let us define by C/3Eiv 
the eigenvalue PDF proportional to 

W \zk-Zjf, Zj = e"'\ -n<ej<TT. (3.40) 

l<j<k<N 

When B = 2 this is the CUE, while the cases B — 1 and Z? = 4 are known in the random matrix literature as the 
COE (circular orthogonal ensemble) and CSE (circular symplectic ensemble) respectively. Similarly, let us define 
by J/?E„ the eigenvalue PDF proportional to 

n 

llx^'{l-Xl)^^ H \Xk-Xjf, 0<Xk<l. 
1=1 l<j<k<n 

In a previous study the two particle distribution functions for the ensembles C/3Ejv in the cases /3 even have 
been expressed as a /9-dimcnsional integral [22], with the latter having the form of an average in the ensemble 
J(4//3)E/3. Generalising the derivation of this identity (see section 3.4) gives 

N m 

In the case /3 = 2 it then follows from (1.21) that for n £ Z>o 

(f\z\-''-'^^'^\l + ziy'+^^{l + tziY) o,{f\z\-''+^''^^'^\l + ziy'{l + {l-t)zif) . (3.42) 

\ / CUE^ \ / CUE^ 

As we have characterised the averages in (3.42) in terms of solutions of the a form of the Pyi equation (1.27), it 
must be that both sides of (3.42) satisfy (1.27) with the same parameters b. 
To check this, we note that according to (3.7) 

UN{t;-N -ri2,n,r)x +772) 

, JV 
= t(t - l)-f log ( it - l)4lbl-i<=2lb]^ie2[b] / TT (,i-^2)/2|;^ ^ ZlV'^+'^Hl + tZlV) ) (3.43) 

at \ \ ^ / cuen/ 

satisfies (1.27) with 

. (N + rj2-Ji N + r,2 + n N-r,2 + n N + r,2+H \ h h h ^ 

b= 2 '^^^ 2 ' 2 ' 2 / V'^MMM)- 

Now the right hand side of (3.42) is obtained from the left hand side of (3.42) by the mappings 

t^ 1-t, 1^ N, rji r]i + r]2, ??2 i-> -m 
which when applied to (3.43) tells us that 

U}f{l-t;-tJ. + ij2,N,rii) 

, N 
= -t{l - t)^ log ({t - l)^<=2[bl^e^[b]-ie2[b] / TT ^(vi-V2)/2^^ _^ Zll^'+^'^l + tZlf) ) (3.44) 
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satisfies (1.27) with b replaced by 

b := (—61, 62, — &4, —bs). 

But we have already commented (recall the paragraph inchiding (3.6)) that —h{l — t) satisfies (1.27) with the 
sign of an odd number of the 6's reversed. Furthermore (1.27) is symmetric in the 6's so we deduce t{t — 1) times 
the logarithmic derivative of the averages in (3.43) and (3.44) satisfy the same equation provided 

te'^lh] - ie2[b] = te^[b] - e'^[h] + ^e2[h], 

which is readily verified. 

Another manifestation of the duahty relation is that for ^ = and n a positive integer, Ejf{t;a,b,n;$,) as 
specified by (1.10) can be written as a /j. x fi determinant. 

Proposition 18. Let T3[yu](i), /j G Z>q, denote the r-function sequence (2.56) with 

T3[l](t) oc 2Fi{-N, N + a + b+l;l + a;t). 

Then 

Eif(t;a,b,iJ.;0) oc T3[n](t) oc t'^'-"''^^^^ det \2Fi(-N - j, N + a + b + 1 + k;l + a;t)] 

I J i,fc=0,...,M-l 

Proof. When = 1 we have 

N 

Eif{t;a,b,l;0) = (Y\{t-xi)) oc Pj^''''' {1 - 2t) oc 2Fi{-N, N + a + b + 1;1 + a;t) (3.45) 

\ p-^ / JUE 

where the first proportionality follows from a standard result in [63] {P^'''-' denotes the Jacobi polynomial of 
degree N). RecaUing Proposition 3, we therefore have 

Ei{t;a,b,l;0)cKT3m;h), b = ( - iV - i(a + 6), -i(6 - a), iV + i(a + 6) + 1, i(a + 6)) . 

Thus 

T-3MW =T-3(t;b), b = (-61,-62,-64,-63) (3.46) 

where b = (61, 62, 63, 63) is specified by (3.9). The operation of reversing the signs of all the 6's and interchanging 
63 and 64 leaves e2[b], e2[b] and the differential equation (1.27) unchanged. It follows from (2.7) that T3(t;b) 
itself is unchanged, so in fact 

T3[n]{t)^r3it;h), 

which is just the equation (3.9). That the case ^ = of (3.9) is singled out by (3.46) follows from the latter being 
a polynomial in t. □ 



3.4 Relationship to generalised hypergeometric functions 

In this section we will show that E'lf{s; a,b, fi) can be identified as an integral representation of a generalised 
hypergeometric function [42] based on Jack polynomials evaluated at a special point. To present this theory 
requires some notation. Let k := (ki, . . . , kn) denote a partition so that Ki > Kj {i < j) and Kj € Z>o. Let 
denote the monomial symmetric function corresponding to the partition k (e.g. if k = (2, 1) then to-k = Z1Z2+Z1Z2), 
and for partitions |«;| = define the dominance partial ordering by the statement that k > ^ if k 7^ /i and 
Xi^=i '^j ^ ^^=1 ^'^'^ each p = 1, . . . , N. Introduce the Jack polynomial Pi^^°' (zi, . . . , zn) =■ Pk^''"^ (z) as the 
unique homogeneous polynomial of degree |k| with the structure 

(the aKfj. axe some coefficients in Q(a)) and which satisfy the orthogonality 

(Pi^/"\Pi^/"V"'oc5.., 
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where 

"1/2 /■1/2 



dxi--- i dxN f{zi,...,ZN)g{zi,...,ZN) TT 

-1/2 J-i/2 l<j<k<N 



Zk Zj\ , Zj . — 6 



We remark that when a = 1 the Jack polynomial coincides with the Schur polynomial. Introduce the generalised 
factorial function 

r„i(«) _ TT - (i - 1)/q + Kj) 

r(«-(j-i)/a) • 

Let 

n + + j))' 

(i,j)eK 

where the notation (i, j) G k refers to the diagram of k, in which each part k, becomes the nodes 1 < i < Kj 

on a square lattice labelled as is conventional for a matrix. The quantity a{i,j) is the so called arm length (the 
number of nodes in row i to the right of column j), while is the leg length (number of nodes in column j 

below row i). Define the renormalised Jack polynomial 

Ci"\z) := ^-plp("\z). (3.47) 

Then the generalised hypergeometric function pFg"^ based on the Jack polynomial (3.47) is specified by the series 

,Ft\au ...,a,-M,..., z) := ^^J^K:^^'"' ^^'^'^ 

(when N = 1 this reduces to the classical definition of p-Fq). 

The relevance of the generalised hypergeometric functions to the present study is that it was shown in [21] 
that 



\Xj - Xk\ 



^ r x^' (1 - xi)^^ {1-txiy--- r dxN x%' (1 - xn)^^ (1 - txNy n 

^ ■^0 •^0 i<j<ks.- 

N 

:=(rr(i-teo-'") 

\ r=i j(2/«)Eiv 

= 2Fi°'Hr,-{N-l) + \i + l;-iN-l) + Xi + \2 + 2;ti,...,tN)\ (3.49) 
a a lti=--.=tN=t 

In the case q = 1, Ai = a, A2 = /j. r = —6, the integrand in (3.49) coincides with the integral in the second 
equality of (1.6). Hence, by Proposition 12 wc have that 

T3{t;h) = 2Fi^\-b,N + a;2N + a + ix;ti,...,tN)\ (3.50) 

where b is specified by (3.2) and T3 is normalised so that T3(0;b) = 1. This, after inserting the known value of 
Pi^' evaluated with all arguments equal in (3.48), can seen to be in agreement with a conjecture of Noumi et 
al. [49] subsequently proved by Taneda [64] using a different argument to that here. 

For future reference we note that for —r =: m, m < N , a non-negative integer, a result of Kaneko [42] gives 
that the average in (3.49) has the alternative generalised hypergeometric function evaluation 

JV 

fj^ / J(2/a)Ejv 

= ^F^^^"^ i-N, -(iV - 1) - a(Ai + 1); -2{N - 1) - a(Ai + A2 + 2); ti, . . . , t„)| 

ltl = '--=tm=t 

0C2Ff'/"'(-jV,-(Ar-l)-Q(Ai + l);Q(A2+m);l-ti,...,l-im)| _ _ _ (3.51) 
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where the final formula follows from an identity in [21]. Note that the role of m and N is interchanged relative 
to (3.49). In fact it is the equality of (3.49) and (3.51) which implies the duality relation (3.41). Also from [42] 
we have that 

JV 

(Y\{t-xir) 0C2Fi^^°'\-N,a{Xi + X2 + m + l) + N -l;a{Xi+my,ti,...,tm)\ ■ (3.52) 

1=1 /J(2/a)Ejv Ul = . ..=t„=t 

In the case a — 1, and with Xi = a, \2 = b, m — fi this coincides with the definition (1.10) of Ej^f{t;a,b, fj,;0). 
Hence from Proposition 13 we have a formula similar to (3.50) valid for n,N£ Z>o 

T3{t;h)oc2F^^\-N,N + a + b + ix;a + iJ,;ti,...,ti,)\ . (3.53) 

Results relating generalised hypergeometric functions to r- functions can also be found in [1]. 



4 Pv scaling limit 

4.1 Scaling limit of the cJUE at the spectrum singularity 

We now turn our attention to the scaled limit (1.26). To anticipate that it is well defined, one notes that the 
integrand defining the average has an interpretation in classical statistical mechanics as a log-gas. Thus it can be 
written in the form of a Boltzmann factor e~'^^ with inverse temperature (3 = 2 and potential energy 

JV 

!7 = -^(^^j+a;ilog|l + ^i| + Mlog|e''"-''/'^'-«i|) - ^ log\zk-z,\. (4.1) 

j=l l<j<k<N 

The classical particle system corresponding to (4.1) has N identical mobile charges interacting via a repulsive 
logarithmic potential, and subject to the discontinuous external potential —{u!2/2)9i {—n < < tt). The mobile 
charges also interact with impurity charges at ^ = tt (of strength uJi) and at ^ = tt — X/N (of strength fi). In 
random matrix language these impurities correspond to spectrum singularities of degenerate eigenvalues. We 
anticipate (1.26) to be well defined because changing variables 9i i— > 9i/N gives a system which has spacing 
between particles of order unity, and X is measured in units of this spacing. 

In the notation of (3.28), the average in (1.26) is denoted ^jv(t;a)i,a;2,A*;C*). It follows from (3.29), (3.34) 
and Proposition 16 that 

t{t - l)^logAN{t;u;i,u>2,n;0 = <^{t) - Git + C2 
where a{t) satisfies (1.27) with b given by (3.31) and Ci, C2 given by (3.33). Since with t = e^^'^ , 

, d d 

'^'-'^dt-''dx 

it follows from the definition of u in (1.26) that 

hm ({ait)-Cit + C2)\ ) =w(X;(a;i,a;2),M;0- 

This limit can be taken directly in the differential equation (1.27). 

Proposition 19. Let Ci, C2 be defined as in (3.33), let b be given by (3.31), and replace h in (1.27) by w+Cii— C2. 
Then with the change of variable t = e''^^^, as iV —> 00 the leading terms of (1.27) are of 0{N^), and give that 

h{t) = u{U; {uJi , W2), m; + + 2CJ1H + ^1/2 (4.2) 

satisfies (I.42) with 

_ _ - itjJ2 . , , „\ 

vi=H, ^2 = -/i, vz=uj, V4 = -uj, Vj:=Vj + —, uj := uJi -\- iijJ2. (4.3) 
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Proof. Direct substitution with the change of variables t = e^-'^l^ and expanding in N for N large shows that 
the leading order term is proportional to N'^ . Equating the coefficient of iV^ to zero gives 



,, X . du „ , wi dw 2 , \2 ^ 

+ ^UJ2 + —)u— + 8iJ,{u)iu;2+ IM2 + —X)— -u - 4wi/iM + (2/ia;2) =0. 
Z dJi. A dX 

Changing variables X i-^ it and comparing the resulting equation with (1.42) after the replacement hy i— »■ 

h + ait + 02 shows that the equations for u(it) and h coincide provided 

1 • 1 2 „ 

ffll = 2*'^2, 02 = ^0J2 + 2Wl/i 

and the equations 

ei[v]=2ioj2, e2[v] = -(m^ +a;? + w|) eslv] = -2ii/ui2, e4[v] = /i^(a;? + a;|) (4.4) 

hold. Direct substitution verifies (4.4) is satisfied by (4.3). □ 
The boundary condition which (4.2) is to satisfy can be predicted in the cases ^ = and ^ = 1 from the 

log-gas interpretation of the average in (1.26). Consider first the case ^ = 0. Then the potential energy function 
(4.1) is being averaged over the whole circle. For large X (after the scaled limit has been taken), to leading order 
the impurity charge originally at ^ = e'^'^"'^^^' on the circle is expected to decouple from the impurity charge at 
z = e"^, meaning that the average in (1.26) will be independent of X. However, for this to happen the potential 
energy must be modified to include the term — 2a;i^log |1 + e^''^^^\ corresponding to the interaction between the 
two impurity charges (see [18] for similar applications of this argument). This leads to the prediction 

2 

u{X;{uii,u!2),fJ-;0) ~ -2uiiiJ. and thus h{t) ~ (4.5) 

X— »oo t— > — zoo 2 2 

In the case ^ = 1, the average (1.26) corresponds to excluding the particles from the interval (tt — X/N,n). To 
leading order one would expect the large X behaviour to bo independent of the parameters wi,W2,/U and thus to 
be given by the uji = uj2 = n = results [69, 17] which implies 

M(X;(a;i,a;2),//;l) ~ -^+0(1) and thus hit) ~ ^ + ^t + 0{l). (4.6) 

X— *00 lb i— > — 2C50 io Z 

4.2 Multi-dimensional integral solutions 

As commented in the Introduction, the differential equation (1.42) is the analogue of (1.27) for the PV system. 
Specifically, one considers [55] a Hamiltonian Hy associated with the Painleve V equation, 

tHv = q{q - 1)V - |(W2 - vi){q - if - 2{vi + V2)q{q - 1) + tq^p + (v3 - Vi){vi ~ Vi){q - 1) 

where the parameters vi, . . . ,V4, are constrained by 

Vl + V2 + V3 + V4 ^ 0. 

It is well known that eliminating p in the Hamilton equations (2.2) gives the general Py equation with 5 = —1/2 
in the usual notation. Our interest is in the further fact that the auxiliary Hamiltonian 

hv = tHv + {v3 — vi){v4 — vi) — vit — 2vi 

satisfies (1.42), which is the analogue of the fact that (2.7) satisfies (1.27). Using the Hamiltonian formalism, we 
have previously constructed determinant /multi-dimensional integral solutions of (1.42) [30]. To make use of these 
results, we first note that with 

hv =a-vit- 2vl (4.7) 
it follows from (1.42) that the quantity cr satisfies 

{ta"f -[a- ta' + 2{a'f + (i^o +v-,+V2 + vz)crf + 4(i/o + a')(v-, + a'){v2 + <r'){v3 + o"') = (4.8) 
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with 

1^0=0, Vl=V2 — Vl, l'2=V3— Wl, l'3=V4 — Wl. (4.9) 

Now, a family of multi-dimensional integral solutions of (4.8) found in [30, Prop. 3.1 together with (1.45)] states 
that 

t|log(r"+"'(e-*S?-^Oj„^J =[/„^(t;a,6) (4.10) 
where U^{t;a,b) satisfies (4.8) with 

vq =0, z/i = —h, z/3 = n + a, 1^3 = n. (4-11) 

The relevance of (4.10) for the present purposes is that in the special case m € Z"*" it follows from (3.42) that 

, N 

lim X-^log(n^^"''''^'|l + «ir'+'''(l+e"''^'^'«0'") =^T^^log(e"''^^=^'OjuE L -n • 

JV— 00 dX \ ' /CUE„ dX ^ /JUE„| ^l-r,2-™ 

(4.12) 

We recognise the left hand side of (4.12) as a rewrite of the left hand side of (1.26) in the case £, = £,* , where 
^* = for m even and (,* =2 for m odd. Consequently 

d , , -xT'^_,xi\ I f.^^ ,a + b + m .h — a^rn^^\ rn . 

Ai=o 
X2=b 

It thus follows from (4.12) and Proposition 19 that 



^^log(e ^'=^ ^>JUE„|,^_ =w(iX;( ^ ,»__),_; 5* j _ _X. (4.13) 



log(e-*S"=i "^•)juE„ = h{t) + 1(6 - a - 2n)t - ^(o + 6 + n) + 1(6 - af (4.14) 
at 4 2 o 

where h{t) satisfies (4.8) with wi = n/2 + (6 — a)/4 and vq, . . . as in (4.11). Use of (4.7) to substitute for h{t) 
in (4.14) reclaims the fact that (4.10) satisfies (4.8) with parameters (4.11). 

4.3 Generalised hypergeometric function expressions 

The generalised hypergeometric function 2^1'™' has the confiuence property [73] 

lim 2F'r\L, 6; c; x/L) = iF[''\h; c; a), (4.15) 
which follows easily from the series definition (3.48). Applying (4.15) to (3.49) with 1 1— >■ — t/r, r — >■ oo gives [20] 

/TTe-*^') =,Fif")(_(n-l)+Ai + l;-(n-l) + Ai+A2 + 2;-fi,...,-u) . (4.16) 

^l-X /j(2/a)B„ \a a / \ti = ---=t„=t 

In the case a = 1, Ai = a, A2 = 6 we can substitute (4.16) in (4.10) to deduce 

t^\ogt''"+"\Fi'\n + a,2n + a + b;-ti,...,-tn)\ =(7^(t;a,6) (4.17) 

where U^{t;a,b) satisfies (4.8) with 

i>o =0, 1^1 = —6, 1^2 = n + a, vs = n. 
Furthermore, applying (1.20) to the left hand side of (4.16) wc deduce [20] 

n 

/ TT + zif+^'e-'^') oc ( - 6'; -!-(n - 1) + a' + 1; ti, . . . I (4.18) 

\^-^ /c(2/a)E„ V a J\ti=-=t„=t 
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5 Applications 



5.1 Eigenvalue distributions in ensembles with unitary symmetry 

Let a sequence of N eigenvalues Xi, . . . , xn have joint distribution pm{xi, . . . ,xm) and let the support of these 
eigenvalues be an interval /. Let lo C /, and consider the multi-dimensional integral 

En{Io;0 '■= (^j J )dxi---(^j J ^dxNPN{xi,. . . ,xn)- (5.1) 

We see immediately that the probability EN,n{Io) of there being exactly n eigenvalues in the interval 7o is related 
to (5.1) by 

Equivalently 

oo 
n=0 

In particular, in an obvious notation it follows from (1.10) and (1-18) that 

(—1)" f)" I 
Ei^„{{s, 1); a, b) = ^-^—Eiis; a, b, 0; ^ (5.3) 

E'J^^{{n-^,n);{uJi,0J2)) = ^^^E^J{<t>;{oJi,0J2),0;d ■ (5-4) 

According to Proposition 13, and the fact that with s = 1, Ej^ in (5.3) equals unity, we have 

Ei^^iis, 1); a, b) = _ exp I - y {ui,{t- a, b, 0; - e'2 [h]t + -62 [b]) ^^^-^ } |^^^ (5.5) 

where U}f is specified as in Proposition 13 with boundary condition (1.38). The result (1.41) with /j, = gives an 
analogous evaluation of (5.4). 

Suppose now that lo = {s,t). Then the quantity 

^^{j^^j )'^^i'"{j~^J '^dXNPN+l{s,Xi,...,XN) (5.6) 

is related to the probability density pN,n{Io,s) of there being an eigenvalue at s as well as n eigenvalues in the 
interval Jo = is,t). Thus we have 

PnAIo,s) = —^—En{Io;s;0\^^^. 

Now for the JUE with lo = (s, 1) the quantity (5.6) is proportional to s''(l — s)'' times Ejf{s; a, b, 2; ^) as specified 
by (1.10). Writing the proportionality in terms of the integral (1.23), it follows from Proposition 13 that 

PAr,„((3,l);a,&) = (iV + l) ^ ' s (1-s) ^— ^ 

JN+i[a,o) n\ 



En{Io;s;0 ~ (N + l) 



gn 



x^exp{-£ [Ui{t; a, 6, 2; - [h]t + \e2 [b]) j^^^ } |^^^ (5.7) 

where tjjq satisfies the differential equation specified as in Proposition 13. The boundary condition which IJ^ 
must satisfy is given by the following result. 

Proposition 20. We have 

U},{t;a,b,2-0 ~ ei[b]-ie2[b]+0(l-t)+do(l -*)''+' [1 + 0(1 -*)] + ••• (5.8) 

t-»i- z 

where 

. ^ 1 rja + b + N + S) r(6 + jV + 3) 

° ^(6+l)(& + 2)r(6 + 3)r(6 + 4) T(a + N) r{N) ^ ' 

and the terms 0{1 — t) are analytic in 1 — t. 
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Proof. The constant term follows immediately from the requirement that the integrand in (5.7) be integrable 
at t = 1. To deduce the explicit form of the leading non analytic term we note that the analogue of (5.2) applied 
to (5.7) gives 

p:w-i,o((s, l);a, b) + (1 - OPAr-i,i((s, 1); o, 6) 

But it follows from the definitions that 

Pw-i,o((s, l);o,&) ~ Pw,i(s)-/ PAr,2(s,i)dt H (5.11) 

s->l- Js 

pjf_i,i((s,l);a,6) ~ / pi,,2{s,t)dt + --- (5.12) 

where plf „ denotes the n-point distribution function in the JUE with N eigenvalues. Hence, making use also of 

the ansatz (5.8), (5.10) reduces to 

PnA^)-^ / P^?,2(s,t)di H 

J s 

~ jV -^^-;(«' b + 2) ^.^^ _ r ^ _ ^) ^ . . . ^ dolN^N r _ ^^,+3 ^ _ ^ . . . \ (5 ^3) 

s— Jn((1,0) I. + J 

For the first term on the left hand side of (5.13) we note from the definition of p^_i(s) as a multiple integral that 

Pjv,i(s) ~ -'V , , s (1-s) , 

s-*l- Jjv(ffl,0) 

which is in agreement with the first term on the right hand side of (5.13). For the term proportional to ^ we note 
from the definition of PN^^i^jt) as a multiple integral that 

s,t->i- JN(a,o) 

which gives 

I\i2{s,t)dt ~ JV(iV-l) -^^-^("'\+^^ l-5)^''+ %^ o^ ^ 
> ) \ ) jj^{a,h) ^ ' (6+l)(6 + 2)(6 + 3) 

Substituting in (5.13) and making use of (1.25) shows this is consistent with the value of do in (5.9). □ 

The probability pjv.n for the JUE with /o = (0, s) can be similarly characterised. In fact by changing variables 

a;; I— > 1 — in (1.10) we see that 

Pw,n((0, s); a, h) = p^,n((s, 1); a)- (5-14) 
Again with 7o = (s, t) we introduce the quantity 

EN{Io;{s,t);0 ■■= (N + 2){N + 1)(^J^-^ )'^*'' " ' ( / / )dxN PN+2is,t,xi, . . . ,xn). (5.15) 



We have 

pN+2[t) n\ dt," le=i 
where pN,n{Io, {s,t)) denotes the probability density of there being an eigenvalue at s and n eigenvalues in 7o = 
{s,t), given there is an eigenvalue at t {pN+2{t) denotes the eigenvalue density at t). For the CUE the eigenvalue 
density is a constant. Clioosiug (s,t) = (tt — tt), (5.15) is proportional to sin^ ^/2 times E'^{<j); (1,0), 2;^) as 
specified by (1.18). Using (1.41) and (1.24) we deduce 

27r\ CUE (2-kX\ 

J^)PN-2,n[^—) 

= -(iV^ - l)sin^ !L_i_^^exp{ - - (^[/^-_,(cot-;(iV,0),2;0 + 4cot-j de}\^^^. (5.16) 
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The boundary condition can be deduced by adapting a strategy analogous to that used in the proof of Proposition 
20. We find 



t^^;-2 (cot ^ ; (AT, 0) , 2; + 4 cot ; 



~ A,;v._4).„?+0(.„.?)+<|zi)l^(.„-| + 0(.„«^)) ,5.17) 

where are term 0(tan^ |) are odd in while all terms 0(tan'' |) are even in 9. 

The eigenvalue density is another statistical property of the eigenvalues accessible from the average (1.10) (for 

the JUE) and (1.18) (for the cJUE). We must choose fi ^ 2 and C = 0, and multiply (1.10) by x"-{l - xf and 
(1.18) by e"'2'^'|l + e''*|2"i. We then obtain expressions proportional to the definition of the density in an ensemble 
of iV + 1 eigenvalues, pjv+i say. Making use of Proposition 13 it therefore follows that 

p^+,(s) = {N + l)ij^±^s%l - sf exp { (f/^(t; a, b, 2; 0) - ei[b]t + ie2[b]) ^^^} (5.18) 

where U'l, is specified as iu Proposition 13 with boundary condition (1.36). An analogous formula can be written 
down for p'§+i{(t>) using (1.41). 

The term in (5.18) involving the exponential function is proportional to Ej^{s; a, b, 2; 0). According to Propo- 
sition 18 and (3.45) 

P^"''''(l-2s) P^"'''+'^(l-2s) 

On the other hand, from the fact that {Pn'^'''\l — 2s)}„=o,i,... are orthogonal with respect to the weight function 
s''(l — s)*", (0 < s < 1), well known direct integration methods give that 



£'^(s;o,6,2;0) oc s ^ 



(5.19) 



Eji{s;a,b,2;0) oc 



P^^^-''\l-2s) ±p'^^''\l^2s) 



(5.20) 



Using Jacobi polynomial identities, (5.20) can be reduced to (5.19). 

The above results apply to the finite N ensembles. Also of interest are scaled iV — > oo limits of these results. 
As mentioned in the introduction, we are restricting attention of such limits to the spectrum singularity. One 
example of the latter type is 

buik/^^N ,. 27r CUE f2nX\ 
p„ (X):=^lim j. 

For its evaluation we read off from (5.16) that 

and note that the scaled limit of (5.15) for the CUE is proportional to the average in (1.26) with 

a;2 = 0, Wi = = 1. (5.21) 

Thus 



2 

bulk 

P 



W-Y^^^|^-P{/ «(^;(M),l;e)f}|^^^ (5.22) 

where according to Proposition 19 the quantity h{t) = u(it; (1, 0), 1; ^) + 2 satisfies (1.42). Substituting the 
parameter values (5.21) in (4.3) it follows from (4.2) and (1.42) that u itself satisfies 

{su"f + {u- su'){u - sm' + 4 - 4(m')^} - 16(u')^ = (5.23) 

with boundary condition 

u{t, (1, 0), 1; -.-y/ + ^(*') - ^457 + ^(*')) (5-24) 
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where all terms 0{t^) are even, while all terms O(t^) are odd. In the case n = 0, p^^^^{X) is the distribution 
of consecutive neighbour spacings in the infinite CUE, scaled so that the mean eigenvalue spacing is unity. The 
formula (5.22) gives 

pr'iX) = ^X^exp ( 1^"" ''^-''-f -\ t), v= (1,-1, 1,-1). (5.25) 

Previous studies have given two alternative Painleve cr- function evaluations of this same quantity. The first, due 
to Jimbo et al. [39], can be written 



Po"'\X) = ^ exp M ll^-p^dt), v = (0,0, 0,0) (5.26) 



while the second, due to the present authors [28], reads 

3ulk.^^ d ( /-^-^ h{-it--v) 



Po'"'^(X) = - — cxpM ^ ^' ' dty v = (0,0, 1,-1) or (1,-1,0,0) (5.27) 
^ 

(the boundary conditions in (5.26) and (5.27) are not required for the present purpose). In each of (5.25)-(5.27) 
the function /!.(t;v) satisfies (1.42) with v as specified. 

The equality between (5.26) and (5.27) implies the identity 

h{s; (1, -1, 0, 0)) = -1 + h{s; (0, 0, 0, 0)) + s^^MMlMl. (5.28) 

h{s; (0,0,0,0)) 

A direct proof of this identity was given in [71] using the formulas of [16] for hv and h'y in (1.42) in terms of a 
particular Py transcendent and its derivative. Analogous to (5.28), the equality between (5.27) and (5.25) implies 
the identity 

h{s; (1, -1, 1, -1)) = -1 + h{s; (0, 0, 1, -1)) + s^gig^Ml^i^. (5.29) 

In fact both (5.28) and (5.29) arc particular examples of an identity between functions satisfying the a-form of 
the Py equation proved in our paper [28]. With a = CT(f;v) specified as a solution of (4.8) with v given by 

V = ( - i(2Ar + a - m), -^(2A^ + o + 3/i), ^{2N + 3a + /x), i(2iV - o + /x)) (5.30) 
the identity states 

a(t;v)| =-(a+l + 2iV)+f+(a(t;v) + t^^)| . (5.31) 

Now, from (4.7), 

h{t; (1, -1, 1, -1)) = a{t; (1, -1, 1, -1)) - t - 2, h{t; (1, -1, 0, 0)) = a{t; (1, -1, 0, 0)) - t - 2, 
/i(t;(0,0,l,-l)) =a(t; (0,0, 1,-1)), /i(i; (0, 0, 0, 0)) = a(t; (0, 0, 0, 0)) 

To deduce (5.28) from (5.31) we set = 2, iV = -1, a = in (5.30), giving v = (1,-1,0,0) on the left hand 
side, then set = 0, AT = 0, a = giving v = (0,0,0,0) on the right hand side of (5.31). The identity (5.28) 
then follows by making use of the second and fourth formulas in (5.32). Similarly, (5.29) follows from (5.31) by 
choosing fj, = 2, N = —2, a = 2, then jj, = 0, N = —1, a = 2, and making use of the first and third identities in 
(5.32). 

There is an analogue of (5.31) for functions satisfying the tr-form (1.27) of the Pyi equation. This can be 
obtained by substituting in the formula 

PAr,n((s, 1); a, b) = ^£^+i,„((s, 1); a, b) (5.33) 
the exact evaluations (5.7) and (5.5), and taking the logarithmic derivative. We find 

Ui{s;a,b,2;0 = -(a + b + 2 + e'2[h]\ - e^[b]| _ )s + (a + 1 + ie2[b]| " Je2[b]| _ ) 

+ (s; a, 6, 0; + s(s - 1) log (i7i^+i(s; a, 6, 0;0- 62 [b] I s + ie2[b]| ) 

= 1/2 - s + Un+1 {s; a, b, 0; + s(s - 1) ^ log (f/^+i (s; a, 6, 0; C) + (iV + 1 + 1 (a + 6)) ' (s - 1/2) + ^ (a' - &')) ■ 

(5.34) 
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An alternative derivation of (5.34) can be given. Following the method first introduced in [71], and also used 
in [30] to derive (5.31), we use (2.26) and the action 

sib = (6i, &2, 64, bs) 

to note that 

siT^^siT^^{bi,b2, 63 + 1, 64 — 1) = (&i, &2, 63, 64 — 2). 
This means that if we temporarily reorder the components in (1.35), 

(i(6-a),-i(a + 6),i(a + 6)+iV,-i(a + 6)-iV-M) 

(as (1-27) is symmetric in the b's, this does not affect the Lf^, which are defined as solutions of (1.27)), then 

Uiis;a,b,2;0 = siTi'siTi'ui+,{s;a,b,0;O (5-35) 

Using (2.3), (2.7), (2.8), Table 1, (2.24) and (2.32), we can use computer algebra to verify that the right hand 
side of (5.35) agrees with the right hand side of (5.34), thus providing an independent derivation of this result. 

5.2 Relationship to the hard edge gap probabihty in the scaled LOE and LSE 

The Laguerre orthogonal ensemble (LOE) refers to the PDF 

1 ^ 

^ncc?e-^'/' n \xk-xj\, XI >0. (5.36) 

1=1 l<j<k<N 

For a — n ~ N — I (n > N), this is realised as the joint eigenvalue distribution of random matrices A, where 
v4 is an n X A'' matrix with independent, identically distributed standard Gaussian random variables. Scaling by 
xi I— » Xi/AN , and taking the limit iV — > 00 gives well defined distributions in the vicinity of a; = (the hard 
edge) [21]. One particular quantity of interest in this scaled limit is Ei'"'^[s; a), the probability that the interval 
(0, s) is free of eigenvalues. This was evaluated in terms of a Painleve V transcendent in [19] and expressed as a 
T-function for a Py system in [31]. The latter result can be written 

^ log£;^'-'^(a;^ (a - l)/2) = hv{x) (5.37) 

where 

, , 5- / \ 1 2 (a — 1) aia — 1) 

av(x) := xnv(x) + -X —- -x -\ ^ — ■ — - 

V / ^ ' 4 2 4 

satisfies (4.8) with 1 1— > 2a; and 

wi = -V2 = ^(a - 1), V3 = -VA = -^(a + 1). (5.38) 
On the other hand, if follows from Proposition 19 that 

u(2ix- (i(a + 1),0), i(a - 1);C) - \{a - l)x + \a{a - 1) 
also satisfies (4.8) with 1 1— > 2a; and parameters (5.38). Hence 

x^ \ogE\^\x'; {a - l)/2) + ^x" = u{2ix; (l(a+ 1),0), l(a - 1);0, (5.39) 

where we use the symbol = to mean that both sides satisfy the same differential equation. 

In the special case (a — l)/2 = m, m G Z>o, if follows from the definition of _Ei'"''*(s; (a — l)/2) that it is 
analytic in s. With this value of (a — l)/2, the right hand side of (5.39) reads u(2ix; {{m+ l)/2, 0), m/2; ^). But 
we know from (4.13) that u with these parameters, analytic in x, is given by an m-dimensional integral, thus 
leading to the following result. 
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Proposition 21. For m G Z>o, 

^^"^''(s^m) = e-"'/*+'"''/e-^''^?'=i"A I . (5.40) 

\ / JUEmlAi=A2=l/2 

Proof. The only point which remains to be checked is that for u in (5.39) given by (4.13), both sides have the 
same small x behaviour, thus determining the boundary condition in their characterisation as the solution of the 
same differential equation. This is the same as showing that both sides of (5.40) have the same small s behaviour. 
Now, with the scaled density of the LOE ensemble (5.36) denoted pa{X), it's easy to see from the definitions that 



But we know that [26] 



^^'""-^'^(s^s') := i((J„(s))' - Ja+lis)Ja-l{s})- 

For small s the term Ja+i(s)/4s dominates, implying the result 

J „2m+l 

To determine the small s behaviour of the right hand side of (5.40), we note that the JUE^, average with 
\i — X-z — 1/2 coincides with an average over the unitary symplectic group USp(m) upon the change of variables 
\j = i(cos0j + 1). Explicitly (see e.g. [30, eq. (1.26)]) 

\ / JUBmlAl=A2=l/2 \ / UeVSp(m) 

Furthermore, it is known from the work of Rains [59] (see also [8]) that 

, > °° / „ /r,\2n rinv 

/ (s/2)TrU\ ^ -Sp [S/2) r^- 

\ /uevsp(m) ^„ 22"n! (2n-l)!! ^ ' 

n=Q 

where f^Z denotes the number of fixed point free involutions of {1,2,..., 2n} constrained so that the length of 
the maximum decreasing subsequence is less than or equal to 2m. Following [1], the small s behaviour of (5.43) 

follows by noting that from the definition, for m < n wc have — 2"(2n — 1)!! (the number of fixed point free 
involutions without constraint) while fm+im = 2™+^ (2m + 1)!! — 1. Hence 



/ (s/2)Tr£/\ r 2,„ [sfZf^"^^''^ 2(m+2)^\ 

).eusp(^) = 1^ - 2-+nm + l)!(2m + l)!! +^^' ^l" 

The derivative of this expression multiplied by has leading order behaviour (5.41), as required. □ 

Note that substituting (5.42) in (5.40) gives the identity 

E'r''((2sf;m) = e-^'/^(e^^ ^)^eusp(™) • (5.44) 

Consider next the Laguerre symplectic ensemble (LSE) which refers to the PDF 

1 ^ 

1=1 l<j<k<N 

The corresponding scaled hard edge gap probability £4"''(s;a) has been shown in [31] to have the t- function 
evaluation 

Et'''{x^;a+1) = ^[fv{x)+f^ix)) (5.45) 

where 

— \ogfv{x) = hv{x) and —\ogfv{x) = hv{x) 
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(recall the meaning of = from (5.39)). Thus 

fv{x) = E\'''''{x^-{a-l)/2) (5.46) 

and the logarithmic derivative of both r-functions in (5.45) satisfy the same differential equation, differing only 
in the boundary condition. From [31] we require 

Analogous to the results (5.40) and (5.44), fv{x) with (a — l)/2 = m, m £ Z>o, can be written as an average 
over the JUE, or equivalently as an average over a classical group. To see this, we first note from Proposition 19 
that 

u(it-, {u)i,Q),n;C) = uiit; (/i,0),wi;O, 
and use this together with (5.39) and (5.46) to deduce 

^ogfv{x) + = u{2ix; (i(o - 1),0), i(o + 1);0- 
But for (o — l)/2 = m and thus (o + l)/4 = (m + l)/2, we see from (4.13) that 

u(2ia;;(i(m-l),0),i(m + l);n = a^:f log(e''"+''^(e-'=^^r=V-A I ). (5.48) 

2 2 ax \ \ / ,7UE„+i I Ai=A2 = -l/2/ 

Because, as will be checked below, the small x behaviour of (5.48) coincides with (5.47), we can read off the 
expression for fyix) as an average over the JUE. 

Proposition 22. For (a — l)/2 = rri, rn G Z>o, 

f*(a;) = e-"'/®+('"+'^^/e-^^^r=V-A I . (5.49) 

\ / JUE^^i I Ai=A2 = -l/2 

Proof. We have to verify the assertion that (5.49) is consistent with (5.47). Again following [1], we proceed in an 

analogous fashion to the proof of Proposition 21. Now, in the case Ai = A2 = —1/2, the JUEm+i average coincides 
with an average over matrices in 0~^{2m + 2) — (2m + 2) x (2m + 2) orthogonal matrices with determinant +1. 
Analogous to (5.42) we have 

\ / JUB^j_i I Ai=A2 = -l/2 \ / UeO+(2m+2) 



But we know from [1] (see also [67]) that 

2 I 

( e ) = exp — ± — + 0{x ^ ) . 

\ /(760±(i) V 2 Z! ^ ') 



ueo±(i) 
Thus 

"^(^)=^"Pl (2m + 2)! +^(^ V' 
which, recalling (a — l)/2 = m, is consistent with (5.47). □ 
We can express £j"'''^(a;^; a + 1) in a form analogous to (5.44). For this purpose we recall that the joint PDF of 
the eigenvalues from the group USp(Z), coincides with the joint PDF of the eigenvalues from the group 0~ (21 + 2) 
(not including the two fixed eigenvalues A = ±1). Using this fact to rewrite the average in (5.44), substituting in 
(5.46), substituting (5.50) in (5.49), and substituting the results in (5.45) shows 

ip-^Vsf/ (a!/2)Tr(7\ / (a;/2)Tri7\ 



Sr'^(x^ 2m + 2) = ie-^'/^ ( (e^^^'^^^) + (e^^/^)^'') 

2 \\ / UeO+(2m.+2) \ /UeO-(2m+2) 

') 

I (/SO(2m+2) 



\ / (/eOf2m+2') 
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5.3 Last passage percolation 

Johansson [40] has introduced the following probabilistic model. Define on each site of the lattice a 
non-negative integer variable w{i,j), chosen independently with the geometric distribution g{q^) (thus Pr(g(g^) = 
fc) = (1 - q^)q^'', fc = 0, 1, . . . ,)• Let (0, 0) wu/wr (M - 1, iV - 1) denote the set of all weakly up / weakly right 
paths (meaning a sequence {(ir,>)}i=i with ii < i2 < ■ ■ ■ < ii and ji < j2 < ■ ■ < ji) on the lattice starting at 
(0, 0) and finishing at (M — 1, iV — 1), moving only upwards or to the right. A quantity of interest is the random 
variable 

Qwu/wr((^ _ 1, TV - 1); g2) := max V w(i,j) (5.52) 

' ' 7re(0,0) wu/wr (M-l,Ar_l) ^ ^ ' 

(i,j)eir 

which with the integer variables w{i,j) regarded as waiting times represents the last passage time of directed 
percolation paths (0, 0) wu/wr (M — 1,N — 1). The relevance to the present study is that Baik and Rains [6, 4] 
have shown 

Pr(G""/'"-((M -1,N- 1); q^) </)=(!- q^)^^ { 11(1 + q\j)^ {1 + U^j)^)^^^ ■ (5-53) 

J = l CUE; 

Using the results of this paper, (5.53) can be expressed both as a 2^"!^' hypergeometric function and as a T-function 
for a Pyi system. 

Proposition 23. We have 

Pr(G'^/^'{{M - 1, TV - 1); g^) < ij (5.54) 

= (l-g^)^^2F«(-iV,-M;Z;ti,...,t,)| 



= (1 - tw^^^i^'' M;M + N-l-t,,...,l-t^) I 

= (l-9^)--r3(9^b)| (5.55) 

' i,= -M-iV-! 

where the second equality is valid for I > N. 

Proof. The average in (5.53) coincides with the CUE average (1.21) provided 

{N,a',b',n,t) = {l,Q,M,N,q''). 

With a, b specified in terms of a', b' in the fine below (1.21), these values imply (o, b) = (—1 — M, M). The 

evaluations now follow from the specification of the right hand side of (5.53) with the above parameters, and the 
results (3.49) and (3.51). For the normalisation in the second equality we have used the general formula [73] 

2F« (r, AT + Ai; 2iV + Ar + A^; ti, . . . , t^)\ = '^^^'1^'''^ . 

The final equality follows from (3.53). □ 
The second equality in (5.54) is well suited to taking the exponential limit 

q = e~^^^'", I = Ls, L —> oo (5.56) 

in which the discrete geometric random variables w{i,j) in (5.52) tend to continuous exponential variables, with 
unit mean, and each random variable again confined to the lattice sites. 

Proposition 24. We have 

g-"/"'((M- 1,A- l);s) := lim Pr('G"'"''""((Af - 1, A - 1); e"'/^) < Ls) 

_r(i + j) 
ri 

j = 

(5.57) 
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where in the final equality it is assumed M > N. Furthermore 



s± logg"'"/--((M - 1, AT - 1); s) = uUt; M-N,0) (5.58) 
where Uji{t;a,0) satisfies (4-8) with 

1/0=0, fi = 0, 1/2 = N + o, z/3 = N. 
Proof. We see from (1.25) that in the limit (5.56) 

Jn{m-n,o) mn tt r(i + i) 
^ ' Jn{m-n,i)^^ nr{M + i + j)' 

while it follows from (4.15) that 

2Fi^\-l,M-M + N-l-ti,...,l-tN)\ iFi^\M;M + N;si,...,SN)\ 

To obtain the final equality in (5.57) we make use of (4.16) and (1.25) in the equality before. The formula (5.58) 
follows from the first equahty in (5.57) and (4.17). □ 

The final equality in (5.57) was first derived by Johansson [40] using a different method. In words this 
equality says that the limiting last passage time cumulative distribution is equal to the probability of there being 
no eigenvalues in the interval (s, oo) of the Laguerre unitary ensemble with a = M — N (the general Laguerre 
weight being x°'e~^). With this identity cstabhshcd, (5.58) is equivalent to a result first derived in [66]. 

Formulas analogous to (5.53) are known for certain variations of the original Johansson model. One such 
variation involves 0, 1 Bernoulli random variables with distribution b{q^) at each site (thus Pr{b{q^) = 0) = 
g^/(l + g^), Pr(6(g^) = 1) = 1/(1 + g^)). Furthermore, the set of weakly up / weakly right paths is replaced 
by the set of weakly up / strictly right paths (to be denoted wu/sr), meaning a sequence {ij.,jr)}r=i with 
ti < 12 < • • • < ii and ii < ^2 < • • • < ji, or the set of strictly up / strictly right paths (to be denoted su/sr), 
meaning a sequence {v, jV)}r=i with ii < i2 < ■ ■ ■ < ii and ji < J2 < • • • < ji- Thus with S^"/'*'' and B^"/^'' 
defined analogously to (5.52), one has [6, 34, 4] 

Pr(B-"/-((M -1,N- lU) <l)= 11(1 + 1/.0*^(1 - (5.59) 

Pr(s="/-((M - 1, AT - 1); q^) <l)={\- q^)'''' I TT(1 + 1/.?0"'^(1 + ^'^O""^) (5-60) 

\ / \ , / CUE; 



where in (5.60) the contours of integration in the CUE; average must be deformed so that the point z = —\ lies 
inside the contour while z = — l/g"^ lies outside. 

Proposition 25. We have 

Pr(B-/-((M - 1, JV - 1); q^) <l) = ^ I^.mn ^^i'^ W -M;l;t,, . . . , ti)\ ^ 

■M-q';h)\ (5.61) 



(1+^2 



1=0, b=N — M-l 



Pr 



(b="/"((M - 1, iV - 1); q2) < /) = (1 - g')*''^2F<''(iV, M; l;ti,..., ti)\^ 



il-qr'rs{q^;h) . (5.62) 

a=0, 6=iV + M-i 



Proof. The averages (5.59) and (5.60) result from the CUE average (1.21) by setting 



{l,0,M,-N,-q^) wu/sr 
{l,0,-M,-N,q'^) su/sr 



(N,a ,b ,n,t) = { ^ 2' _ ;_ (5.63) 
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With a',&' so specified, it follows from the equations in the line below (1.21) that (a, 6) equals {—I — M,M) and 
{—l + M, —M) respectively. The 2^1^^ evaluations now follow from the specification of the left hand side of (5.53) 
with the above parameters, and the results (3.49) and (3.51), and the r-function evaluations in turn follow from 
(3.53). □ 
We remark that for l> N it follows from the definitions that 

Pr(^B""/""((M -1,N - 1); q^) <l) = Pr(^B""/"'^((M - 1, iV - l)-q^) < /) = 1. 
These equalities iov I = N can be read off from the above 2^1^' evaluations by noting the general formula [73] 

m 

2F«(a,6; a;U, ■ ■ ■ ,tm) = iF^^\h;tu ■ ■ ■ ,tm) = ^(l-^,•)-^ 

i=i 

Underlying the Johansson probabilistic model and its generalisations is the Knuth correspondence between 
integer matrices, generalised permutations and semi-standard tableaux (see [40, §2.1]). In fact the random variable 
(5.52) in the Johansson model is equivalent to the random variable specifying the length of the longest increasing 
subsequence of random generalised permutations with a certain probability measure, or alternatively the length 
of the first row of a certain class of random Young diagrams. The notion of measures on Young diagrams also 
occurs in the work of Borodin and Olshanski on representations of the infinite symmetric group (see e.g. [13, 15]). 
In particular, for random Young diagrams distributed according to a so called ^:-measure, the probability that the 
first row length Y{z; ^) is less than or equal to / is shown to be given by 

Pr(y(«;0 <l)= (1 11(1 + ^f^z.ni^-^^lz.y) . (5.64) 

\ / \ , / CUE? 

3=1 

Note that with z = N', ^ = q^, this coincides with the case M = N oi (5.53). Comparing (5.64) with (1.18), we 
see that in the case z real 

Pr(y(z; e"*) < z) = (1 - (-it; (0, 0), 0). (5.65) 

Hence we have the Pyi cr-function evaluation (3.25). Complimentary to this result, we draw attention to the recent 
work of Borodin [10] (see [11] for still more recent developments), who through a newly developed theory of discrete 
integral operators and discrete Riemann-Hilbert problems, has derived a difference equation for Pr^y(«; e~*) < 
in the discrete variable I involving two auxiliary quantities which satisfy a particular case of the discrete Py system 
(2.71), (2.72). 



5.4 Symmetrised last passage percolation and distribution of the largest 
eigenvalue in the finite LOE and LSE 

Consider the original Johansson probabilistic model introduced in the previous section. Choose the weights 
w{i,j) for i < j as independent geometric random variables with distribution g{q^), and impose the symmetry 
that w{i,j) — w{j,i) for i > j. Choose the weights w{i,i) on the diagonal independently with distribution g{q), 
and consider the random variable (5.52) with M = N. By the symmetry constraint the set of paths in (5.52) can 
be restricted to the triangular region i < j, so we denote (5.52) in this case by G^"/^'ie(-^; 
Baik and Rains [6] have shown that 

Pr (Cr/nlJ^; 1') < = - 1^"'^^ - iT'^'-^'^'idet [(1 + U){1 + lUf])^^^^^. (5.66) 

The significance of (5.66) from the viewpoint of the present study is that it can be written in the form (1.7). To 
see this we recall that for I even all eigenvalues of 0~^{l) come in complex conjugate pairs e**^, while for I odd 
one eigenvalue of 0~^(l) is equal to +1 and the rest come in complex conjugate pairs. Hence 

[1/2] 

det [(l + [/)(l + g?7)'^] = xi,N l[{2 + 2 cos e + q^ + 2q cose j)'^ 

3=1 
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where xi,n = 1, I even, xi,n = 2(1 + q)^ , I odd. Introducing the variables 



A.- = ^(cos^. + l) (5.67) 



this reads 

[1/2] 



det [(1 + f/)(l + qUf] = 2'^™%^^'^\,^ n - if + ^0 



.^1 4« 

Furthermore, in terms of the variables (5.67), the PDF of the eigenvalues {e**^^ }j=i^...^[(/2j for matrices U € 0'^{l) 
is of the form (1.1) with N replaced by AT*, w{x) = a;"(l — a;)** (Jacobi weight) and parameters 

{N*,a,b) = {[1/2], -1/2, {-iy-'/2). 

Thus we can rewrite (5.66) to read 

[1/2] 

A.ff (l-# + A,)") 

where 



Pr(G-iTJAr; ,^) < l) = a,.(,)( n A. (^(1 - + A,)^)^^^ I ^ (5.68) 



C(.^(9) := 2^(^+^^>'/^lg^['/^lx<,Jv^(l - qf{l - qY^"-'^^' (5.69) 
The result (3.52) can be used to express (5.68) as a generalised hypergeometric function. 

Proposition 26. With Ci^niq) specified by (5.69) and Jn{a,b) specified by (1.25) 

p /^wu/w. .^.2. < A _ ^, , .„^ ^^[V2](l/2+jV,(-l)'-V2) 

X2F«(-[«/2],A^+[«/2] + (l + (-l)'"')/2;iV + l/2;ti,...,t^)| (5.70) 

valid for [1/2] > N. 

The expression (5.70) is well suited to taking the exponential limit (5.56). 
Proposition 27. We have 

9Z-I::u{N;s) := £m^Pr(G:,r/nIie(iV;e-^/^) < Ls) 

= n i;^±lI.-(-+i)/^e-~=/^ oi^i'^'(_; iV + 1/2; U,...,tr,) 



.^^r(2i + i)" -1 v-,^- . 



(5.71) 



Proof. Suppose for definiteness that I is even. In the limit (5.56) we can check from the definition (5.69) and 
the evaluation (1.25), making use of Stirling's formula, that 



J[i/2](i/2 + iv,(-i)' V2) JX- r(j + 1) Ar(JV+l)/2 -Ns/4 



^''^^'^ j,/2i(-i/2,(-i)'-V2) lAiWT 



s 



J[i/2](-i/2, (-i)'-V2) r(2i + 1) 

Also, analogous to the result (4.15), we can see from the series definition (3.48) that in the limit (5.56) 

^Fi'H-[l/2],N + [l/2] + {l + i-ly-')/2■,N + l/2■t^,...,t^)\ 

~ o-F'i'^_;iV+l/2;ti,...,tjv)| , . (5.72) 

□ 

In [23] the generalised hypergeometric function 

oFi^\_;N + n;ti,...,tN)\ 

I tl = --- = tj^=t 
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for fj, = and n = 2 has been evaluated as an iV-dimensional determinant. Generalising the working therein gives 
that for general fj, 



oF«(_;iV + M;ti,...,t^)| = 11 ^wtlt/^ (7) '''' ^et (.)] (5.73) 

Ul = .-.=tjV=s2/4 -M: i (j + 1) \S/ L i j,k=0,...,N-l 



j=0 

Furthermore, we know from [30] that 



n ro-+i) (7!) '''' K^-'=+''(^)J,,.=o = ^"Pio ^ — ^'-''^ 



where 



v{t; N,n) = - (auv (t) + m(m + N)/2^ . (5.75) 
The frinction ami (t) satisfies the differential equation 

{taui'f - viV2{aui'f + o-ni'(4crni' - 1)(<Tiii/ - taur) - ^{vi - V2f = 0, (5.76) 

which is the Jimbo-Miwa-Okamoto cr-form of the Pm' equation (a simple change of variables relates Piii to 
Piii' [56]) with parameters 

{vi,V2) = in + N,-n + N) (5.77) 

and subject to the boundary conditions 



t^oo 4 2 V 4 2 
The boundary condition (5.78) does not distinguish solutions with \x 1— > — /i (note from (5.77) that under this 
mapping v\ <-» U2, and this latter transformation leaves (5.76) unchanged). However, it follows from (5.73)-(5.75) 
that 



Collecting together the above results allows the evaluation (5.71) of ; s) to be further developed. 



This boundary condition docs distinguish solutions with [i 

tion (5.71) of 
Proposition 28. We have 



j,k=0,...,N-l 



: exp (- r ( - auv{t)\ + i/4 - Nt'/'/2 + {N' - l/2)/4) ^ ] (5.80) 



In the Okamoto r-function theory of the Pni' system [56], the transcendent —aiiii{t)/t is an auxiliary Hamil- 
tonian, being equal to the original Hamiltonian plus a function of t. Changing the function of t doesn't alter this 
property, so we have that 

9ZrI:;i.{N; s) = TuA{s/4f)\.,=^+,;, (5.81) 

V2 = N-l/2 

where 

4logr„p(t)l =-ai„,(*)L^=„+,/, +V4-iVt^/V2 + (iV'-l/2)/4. (5.82) 

' v2 = N-l/2 ' V2=N-l/2 

The symmetrised Johansson model has been generalised by Baik and Rains [6, 7] to include a parameter a in 
the geometric distribution of the waiting times w{i, i) of the diagonal sites. Thus these waiting times are chosen to 
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have distribution g{aq), with the case a = 1 corresponding to the original symmetrised Johansson model. With 
the random variable (5.52) now denoted G^"£^[g(Ar; a, g^), one then has [4] 

P<G:i/:;jAr;a,g^) < z) = 1(1 - aq^l 

X (( det [(1 + aUKl + ,Uf] )^^^^^^^ + ( det [(1 + af/)(l + <lUr])^^^_J (5.83) 

The analogue of Propositions 26, 27 and 28 can be obtained for the case a = 0. 
Proposition 29. Let xIn — ^' ' ^'ven and xIn = (1 + i ^ odd, and put 

CUl) = 2'""/''9""^''x^.i(l - 5^)-(--^)/^ 
Let XiN = (1 ~ 9^)^; ^ even and Xin = (1 ~ 9)''; ' '^'^'^i '^'^'^ P"* 

C^,(g) = 2^-[('-i)/^lg-[('-^)/^'xz>^(l - g^)'^("-^>/^ 

We have 

p /^wu/w. . 2. ^,\_c+ .^ /[V2l(^l/2 + iV,(^l)'-V2) 

X2Fi(^^(-[V2],iV + [V2] + ((-l)'-^-l)/2;iV-l/2;ti,...,t;v)| 

, ^- ., /[0-i)/2](l/2 + ^.(-l)V2) 
+ ^'--^«^ J,/2, (1/2, (-1)72) 

X2J'i''^(-[(«-l)/2],iV+[(Z-l)/2] + (l + (-l)')/2;iV + l/2;ti,...,t^)| . (5.84) 

valid for [Z/2] > A^. 

Proof. The rewrite of the average over 0'^{l) in (5.83) is done in the same way as in going from (5.66) to (5.70). 
To rewrite the average over 0~ {I) in (5.83) we must first recall that for I even there are eigenvalues ±1, while the 
remaining eigenvalues come in complex conjugate pairs e**^^, and for I odd there is an eigenvalue —1, with the 

remaining eigenvalues also coming in complex conjugate pairs e*'*^ . Furthermore, in terms of the variables (5.67), 
the PDF of the eigenvalues {6=*='®^ }j=i,...,[(i_i)/2j is of the form (1.1) with N replaced by N* , w{x) = x'^{l — x)^ 
and parameters 

(iV*,a,6) = ([(/-l)/2],l/2,(-l)72). 
Using these facts, we proceed as in the rewrite of the 0+ {I) average. □ 

Proposition 30. We have 

9ZiZ.{N-Q,s) := \i^Jr(GltZSN;Q,e-'/'^) < Ls) 



\ { n ^^^f'^^''"e-^'l\F^^^ (_; AT + 1/2; ^„ . . . , t.) 



tl = ---=tjV=s2/64 



JV 



-TT r(j) iV(iV-l)/2 -iVs/4 ./V_1/2-f, 

+ ilr(2j-i)' ^ ^-'^ ^/^'^^'■■•'*^^L,=...=t„=.2/64 



ijexpf- r (-a„p(t)| +V4-iVti/V2 + (iV^-l/2)/4)^') 

+ expf-/°° (_^„j,(t)| +t/4-iVtVV2 + (iV2-l/2)/4)^')| (5.8 
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Proof. These formulas follow from Proposition 29 by following the working which led from Proposition 26 to 
Propositions 27 and 28. The two terms in (5.85) correspond to the two terms in (5.84) but in the reverse order. 

□ 

Baik [5] has recently evaluated the scaled limit of (5.83) for general a (the variable a is also involved in the 
scaling). The evaluation is given in terms of the solution of a Riemann-Hilbert problem. It is pointed out in [5] 
that this has consequence with regard to the distribution of the largest eigenvalue in an interpolating Laguerre 
ensemble, which passes continuously from the LOEjv to the LSE;v/2- This follows from the identity [6] (see also 
[32]) 

9ZiZ.{N:A,s) := £m^Pr(G^,r/nIleWe^^'^e-'''^) < Ls) 

Jo Jo 

where xt = 1 if T is true and xt = otherwise. The case ^ = corresponds to q = 1. Now, with A = the 
integral in (5.86) is precisely the gap probability Ei{(s,oo); N;0) (no eigenvalues in the interval {s,oo)) for the 
LOEjv with parameter a = 0. Equating (5.86) in this case with (5.81) then gives the r-function evaluation 

E]'°^iis,^),N-Q) = ruv{{s/4')\^,^.+,/,- (5-87) 

V2 = N-l/2 

Taking the A — > — oo limit in (5.86) corresponds to a = 0. We see from (5.86) that for AT even 

, lim 9::^IZeiN; A, s) = £^^""{{3, 00); N/2; 0), 

where -£4 ^^((s, 00); n; a) denotes the probability that there are no eigenvalues in the LSE„ with weight x"e^^. 
Equating (5.86) in this case with (5.81) then gives 

S4''^^((s,(X));iV/2;0) = ^ (riii'((s/4)^)l + r„p((s/4)^) I ) . (5.88) 

' „2=iV-l/2 't,2 = -N + l/2^ 

As the differential equation (5.76) is unchanged by the interchange vi •*-» V2 the Hamiltonians for both r-functions 
in (5.88) satisfy the same differential equation. It is interesting to compare (5.88) with the structural formula [27] 

E^s, 00); N/2., 0) = 1 (^-((., 00); N; 0) + ^S^;^) • (5-89) 
Doing this tells us that 

E^'^'^iis, 00); AT; 0) = mv{{s/4f)\ mv{{s/^f)\ . (5.90) 

V2 = N-l/2 „2 = -N + l/2 

On the other hand we know from [65] that 

E^''^iis,^y,N;Q)=e^p(- H a{t)\ y) = rv(.)| , (5.91) 

where a{t) satisfies (4.8), and so 

Similar r-function identities to (5.92) have occured in the works [25, 70]. The Painleve transcendent evaluations 
of E]:°^{{s,ooy,N;0) and £;i'^'=((s, 00); iV/2; 0) given in [5] differ from (5.87) and (5.88), involving instead of 
T-functions for the PHI' system, the transformed PV transcendent found in [65] to be simply related to the 
derivative of <j(t)\ 

V ^1 1/0 = 1/1=0 
<^2 = ''3 = '^ 
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5.5 Diagonal-diagonal correlation in the 2d Ising model 

The square lattice Ising model has one of two possible states, Oij = ±1, on each site of the two-dimensional 
square lattice (see e.g. [9]). The infinite square lattice of states is achieved by considering a sequence of finite 
lattices, of dimension (2iV + 1) x (2A'^ + 1) say, centred about the origin. The joint probability density function 
for a particular configuration of states in the finite system is given by 

^ N N-1 N N-1 

— exp ^ ^ Gijai+^j-^-K-i ^ ^ (TijO-ij+i] , 

where Z2n+i is the normalisation. Thus there is a coupling between nearest neighbours in the horizontal and 
vertical direction. In the limit iV — » oo, an unpublished result of Onsager (see [44]) gives that the diagonal 
spin-spin correlation has the Toeplitz form 

(o'oofnn) = det[ai_3]ij=i,...,„ (5.93) 

where 

''^ = h r^^^^""1 V- aAy^T ^'' fc:=«inh2ifisinh2if.. (5.94) 
Regarding the Fourier integral (5.94) as a contour integral, deforming the contour of integration and then writing 
the Toeplitz determinant as an average of the CUE we see that 

i„ „ N_j (nr=i^z^'ii+^r'/'(i + (i/fc')-o'^')cuE„, i/fc^<i 

^ {[ll^l^l ll + '^il ' (l + KZl)' )CUE„, fe <1 



We can now read off from (3.7) the following results. 
Proposition 31. Let s = We have 

+ n)s + i + s(s - 1)^ log(cToo<T„„) = U^{s; 1/2 - n, 1/2, -1/2) (5.96) 

where satisfies the Pyi equation in a-form (1.27) with parameters 

/I/ n 1 , n\ 

^=(2("-^)'2'2(" + ^)'-2)- 

Let t = . We have 

2-1 J 

-^t - g + t(t - 1)^ log(<TOO(Tnn) = Ui{t; -1/2 - 71, 1/2, -1/2) (5.97) 



where ?7„ satisfies the Pyi equation in a-form (1.27) with parameters 

f n n — \ n n + 1 
b = - 



,2' 2 '2' 2 
It is straightforward to check from (5.97) and (1.27) that 

a{t) := t{t - 1)^ log{(T00O"nn) - ^ 

satisfies the differential equation 

('<• - »$)^ = "■((' - - - *w(<' — i) h'i - 

This is a result due to Jimbo and Miwa [37]. Our point therefore is not a new characterisation of {aooffnn), but 
rather the fact that the known characterisation fits our development of the Okamoto r function theory of Pvi . 
We note too that a recent result of Borodin [10] can be used to give a recurrence for {aooann) as a function of n. 
This recurrence applies to all Toeplitz determinants 

ql^'^ := (1 - 0" det[5i_j]ij=i,...,„ 

>e-'f«(l - v^|e'*)"(l - v^e""*)^'. (5.98) 



-/ 
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Comparison with (5.93) and (5.94) shows 

/ \ ■'■ (-1/2,1/2,1/fc^) 

Recent work of Adler and van Moerbeke [2] gives a different recurrence relation for (5.5) to that of Borodin, 
and thus a further recurrence for {(tooo"„„). Explicit computation of {(7ooO"„„) as a power series in , for which 
both the differential and difference equation characterisations are well suited, is a crucial step in recent very high 
precision numerical studies of the zero field susceptibility [47, 48, 57, 58]. 
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